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Universitetskii pr., 13, Moscow, 119992, RUSSIA 



Abstract 

In this paper we review the field-theoretical approach. In this framework pertur- 
bations in general relativity as well as in an arbittary D-dimensional metric theory 
are described and studied. A background, on which the perturbations propagate, is a 
solution (arbitrary) of the theory. Lagrangian for perturbations is defined, and field 
equations for perturbations are derived from the variational principle. These equa- 
tions are exact, equivalent to the equations in the standard formulation and have a 
form that permits an easy and natural expansion to an arbitrary order Being covari- 
ant, the field-theoretical description is also invariant under gauge (inner) transforma- 
tions, which can be presented both in exact and approximate forms. Following the 
usual field-theoretical prescriptions, conserved quantities for perturbations are con- 
structed. Conserved currents are expressed through divergences of superpotentials — 
antisymmetric tensor densities. This form allows to relate a necessity to consider local 
properties of perturbations with a theoretical representation of the quasi-local nature 
of conserved quantities in metric theories. Properties of the conserved quantities un- 
der gauge transformations are established and analyzed, this allows to describe the 
well known non-localization problem in explicit mathematical expressions and oper- 
ate with them. Applications of the formalism in general relativity for studying 1) the 
falloff at spatial infinity in asymptotically flat spacetimes, 2) linear perturbations on 
Friedmann-Robertson- Walker backgrounds, 3) a closed Friedmann world and 4) black 
holes presented as gravitationally-field configurations in a Minkowski space, are re- 
viewed. Possible applications of the formalism in cosmology and astrophysics are 
also discussed. Generalized formulae for an arbitrary metric /^-dimensional theory 
are tested to calculate the mass of a Schwarzschild-anti-de Sitter black hole in the 
Einstein-Gauss-Bonnet gravity. 
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1 Introduction and preliminaries 

1.1 Perturbations in gravitational theories, cosmology and relativistic astro- 
physics 

Much of research in general relativity (GR) is frequently carried out under the assump- 
tion that perturbations of different kinds propagate in a given (fixed) background spacetime 
(exact solution to the Einstein equations) [1] - [4]. A majority of cosmological and astro- 
physical problems are also studied in the framework of a perturbation approach. It is quite 
impossible to give a more or less full bibliography on this topic. Nevetherless, to stress the 
importance of such studies, we shall outline shortly some of the related directions. 

Cosmological perturbations on Freidmann-Robertson-Walker (FRW) backgrounds have 
been considered beginning from the famous work by Lifshitz [5]. The Lifshitz principles 
were developed by many authors in a variety of approaches (see, e.g., the popular review 
[6], and a recent review [7]). As well known examples, one can note the following. The 
effect of amplification of gravitational waves in an isotropic world was discovered by Gr- 
ishchuk [8] ; the gauge invariant theory of perturbations was formulated by Lukash [9] and 
Bardeen [10]. Lately, non-trivial perturbations in FRW worlds are being considered more 
frequently. For example, note the recent papers [11, 12] where, using the quasi-isotropic 
expansions, the authors describe non-decreasing modes of adiabatic and isocurvature scalar 
perturbations, and gravitational waves close to cosmological singularity. Deviations of a 
space-time metric from the homogeneous isotropic background become large; while lo- 
cally measurable quantities, like Riemann tensor components, are still close to their FRW 
values. An approach, where integrals of FRW models are presented in a time independent 
form both for "vacuum" and for "usual" matter, has been developed [13]. In the two last 
decades, cosmological perturbations are considered not only in the linear approximation, 
but including the second order also (see the recent review [14]). 

The evolution of quantized fields, including gravitational and electromagnetic fields, 
is intensively studied on curved backgrounds of various general classes, such as globally 
hyperboUc, static with symmetry groups, etc. In this framework, exact solutions, Uke FRW, 
anti-de Sitter (AdS) and Bianchi of different types, are also exploited as backgrounds. As an 
example, note the theory which is being developed by Grishchuk (see the reviews [15, 16] 
and references there in). Relic gravitational waves and primordial density perturbations are 
generated by strong variable gravitational field of the early FRW universe. The generating 
mechanism is the parametric amplification of the zero-point quantum oscillations. These 
generated fields have specific statistical properties of squeezed vacuum quantum states. 
Cosmological perturbations at the early inflation stage [17] also continue to be studied 
(see, e.g., recent papers [18] - [20] and, e.g., reviews [21] - [23] and references there in). 
Separately, the great interest to AdS spaces has been initiated by the discovery of the present 
(not at the earlier inflation stage) accelerated cosmological expansion (see review [24]). To 
explain the acceleration new cosmological solutions are searched [25]. On backgrounds of 
such solutions perturbations are also to be considered. 

Several black hole solutions [26], which could represent the neighborhoods of relativis- 
tic astrophysics objects, also play a role of backgrounds for evolution of different kinds 
of perturbations. Amphfication and dispersion of metric perturbations (including gravita- 
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tional waves), electromagnetic field and massless and mass scalar fields are studied on these 
backgrounds (see, e.g., resent works [27] - [33]). 

The rapid development of the detecting technique stimulates a development of the 
gravitational-wave physics (see the recent detailed review [34] and references there in). 
Thus the theoretical study of propagation and interaction of gravitational waves becomes 
especially important. In the works by Alekseev (see [35, 36] and references there in), the 
monodromy transform approach was developed for constructing exact solutions of the Ein- 
stein equations with spacetime isometrics. This method has been applied for exact solving 
the characteristic initial value problem of the collision and subsequent nonlinear interaction 
of plane gravitational or gravitational and electromagnetic waves with distinct wave fronts 
in a Minkowski space [37] - [39]. 

Currently, following an extraordinary interest with brane models, D-dimensional metric 
theories of gravitation have been examined more and more intensively. Fundamental works 
on brane worlds have appeared two decades ago (as we know the works by Rubakov and 
Shaposhnikov [40] and by Akama [41] are the first in this direction). Later, especially 
after the works by Randall and Sundrum [42, 43], the interest to these models has risen 
significantly (see, e.g., a review [44]). Perturbations, including gravitational waves, in the 
framework of -dimensional metric theories and brane models are also studied intensively 
(see, e.g., [45] - [51] and references there in). 

The list of directions, where perturbations on curved or flat backgrounds are studied, 
could be continued. Impressing results were obtained in this direction. However, we would 
like to note the problems associated with methods of such investigations, rather than the 
results. As a rule, these methods are restricted since they are constructed for examination 
of particular tasks only. Thus: 

• Although the modern cosmic experimental and observable data require more detailed 
theoretical results, frequently studies are carried out in the linear approximation only, 
without taking into account the "back reaction". 

• Applying particular methods, one generally uses many additional assumptions. Thus, 
it is not clear: what results are more general and what results change under a change 
of these assumptions. 

• It is difficult to understand: could an approach developed for a one concrete back- 
ground be applied to other backgrounds. Frequently only simplified backgrounds are 
used; etc. 

• In each particular case, considering perturbations, one needs to fix gauge freedoms. 
A concrete fixation is connected with a concrete mapping of a perturbed spacetime 
onto a background spacetime. It turns out, that it is not so simple to understand what 
gauge is "better", or how to find a connection between different gauges, etc. 

1.2 Conservation laws and their properties 

Very important characteristics for studying perturbations are such quantities as energy- 
momentum, angular momentum, their densities, fluxes, etc. However, as is well known, the 
definition of energy and other conserved quantities in GR has principal problems, which are 
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well described in many textbooks (see, e.g., [3]). We repeat some related issues to stress 
importance of these notions in a theoretical development of gravitational theory also. 

It is useful to reconsider works by Einstein who paid special attention to conservation 
laws for "energy components" t/j," of the gravitational field. From the early stage, when 
the theory was not yet presented in a satisfactory form (its equations were R^i, = kT^^), 
he examined the conservation law ^^(t^u'^ + T^'^) = [52, 53]. Then, the final form 
R^jj = K [r^y — ^QixyTct'^^ of the GR equations was given in the work [54]. Einstein 

himself explains this change [54] by saying that only the additional term — ^g^vTa^ leads to 
a situation, when energy complexes both for the gravitational field, t^'^ , and for matter, T^^, 
enter the field equations in the same manner. Thus, historically the analysis of conserved 
quantities and conservation laws was crucial for constructing GR. 

In the work [55], Einstein finally suggested the canonical energy-momentum complex 
for the gravitational field t^^. Later it was called as the Einstein pseudotensor. As an appli- 
cation, Einstein used t^i'^ to outhne gravitational waves [56] - [58]. From the beginning he 
stressed that tf/ is a tensor only under linear transformations. Under the general coordinate 
transformations, t^/ can change and even be equal to zero. Einstein interpreted this as a 
"non-localization" of gravitational energy, which is a special property of the gravitational 
field, and not a defect of the theory. This was a reason for numerous criticism and discus- 
sions. Protecting his theory, Einstein himself was the first who gave physically reasonable 
arguments (see, e.g., [59]). 

The criticism far from kilhng the theory, was a reason for the intensive study of prop- 
erties of gravitational field in GR and its further development. Following Einstein, in next 
decades a great number of methods were suggested for defining conserved quantities in GR. 
As a result of these efforts, important theoretical tests, which restrict an ambiguity in the 
definition of conserved quantities, were elaborated. Thus, mathematical expressions have 
to give acceptable quantities for black hole masses, for angular momentum in the Kerr so- 
lution, for fluxes of energy and momentum in the Bondi solution; also, a positive energy 
density for weak gravitational waves. However, due to nontrivial pecuUar properties of 
conserved quantities in GR, up to now: 

• Frequently, it is very difficult to find a connection between different definitions; 
sometimes definitions even contradict one another; sometimes definitions (especially 
earlier) are not-covariant; etc. 

• Sometimes definitions are not connected with perturbations. 

Now, let us present the modern point of view on the non-localization problem. Con- 
sidering the physical foundation of the theory it is clear that the non-localization is directly 
connected with the equivalence principle (see, e.g., [3]). On the other hand, the situation 
can be also explained by the fact that GR is a geometric theory where spacetime, in which 
all physical fields propagate, itself is a dynamical object. (Of course, the non-localization 
problem is related to all the mettic theories of gravity, not only to GR.) Due to these objec- 
tive reasons, sometimes the problem of conserved quantities in metric theories is presented 
as ill-defined. However, only the fact of non-locafization cannot imply that these notions 
are meaningless. Without a doubt gravitational interaction contributes to the total energy- 
momentum of gravitating system [3]. Indeed, describing a binary star system one needs 
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to consider gravitational energy as a binding energy; considering gravitational waves in an 
empty domain of space one finds a positive energy of this domain as a whole; etc. All of 
these are related to non-Iocal characteristics. This conclusion is supported by the mathe- 
matical content of GR. Szabados [60] clearly expresses it as follows. "... the Christoffel 
symbols are not tensorial, but they do have geometric, and hence physical content, namely 
the linear connection. Indeed, the connection is a non-local geometric object, connecting 
the fibers of vector bundle over different points of the base manifold. Hence any expression 
of the connection coefficients, in particular the gravitational energy-momentum or angu- 
lar momentum, must also be non-local. In fact, although the connection coefficients at a 
given point can be taken zero by appropriate coordinate/gauge transformation, they cannot 
be transformed to zero on an open domain unless the connection is flat." 

Thus the non-localization is natural and inevitable. However, up to now there is no 
simple and clear description of it. Therefore: 

• It is important to give mathematical expressions, which constructively present the 
non-localization of conserved quantities in metric theories. 

The non-localization has to be connected with gauge properties in the description of per- 
turbations. Thus, a solution of this problem has to give a possibility to find reasonable 
assumptions for a gauge fixation, which allow to describe certain quantities as localized 
ones. 

It is natural that due to the aforementioned pecuhar properties of gravitational field 
much attention has been paid just to non-local characteristics. Thus the total energy- 
momentum and angular momentum of a gravitating system in whole spacetime were studied 
intensively. Such quantities are frequently called as global ones. In this context asymptot- 
ically flat spacetimes are considered in details (see, for example, earher reviews [61, 62], 
also recent papers and reviews [63] - [67] and numerous references there in). One of the 
great achievements was the proof of the positivity of the total energy for an isolated system 
[70] - [73] (see also the review [74]). The global conserved quantities for asymptotically 
curved backgrounds (like AdS space and some others) are also studied intensively (see, e.g., 
[68] - [78]). 

The aforementioned development has initiated a more intensive examination of the en- 
ergy problem in GR. Conserved quantities became to be associated with finite spacetime 
domains. Such quantities are called as quasilocal ones and can give a more detailed infor- 
mation than the global quantities. In the last two-three decades the quasilocal approach has 
became very popular. It is not our goal to present it here, moreover, recently a nice review 
by Szabados [60] has appeared. Nevetherless, below we shortly outUne some of important 
quasilocal methods. 

The Brown and York approach [79] is based on the generalized Hamilton- Jacobi analy- 
sis. It considers a spatially restricted gravitating system on 3-dimensional spacelike section 
E. A history of the boundary is a 3-dimensional timelike surface S (cylinder). It is as- 
sumed that a 3-metric jij on S is fixed and plays a role of a time interval in the usual 
non-relativistic mechanics, which defines initial and final configurations. They define the 
energy-momentum tensor t^^ on S, as a functional derivative of an action with respect to 
7ij. An intersection of E with 5* is a 2-sphere B which is just a spatial boundary of the 
system. Normal and tangential projections of r^^ onto B give surface densities of energy. 
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momentum and space tensions on B which are quasi-local expressions. It is crucial to deter- 
mine a reference flat space, which is uniquely defined by the isometric embedding B (with 
a positive inner curvature) into a flat space. The Brown- York method received a significant 
development in the works by Brown, Lau and York [80] - [85] and in works of other authors 
(see review [60]). The recent work [86] could be considered as a mathematical textbook on 
this approach. 

At the first stages of constructing the Hamiltonian dynamics of GR by Arnowitt, Deser 
and Misner [87] (ADM), surface integrals were neglected a priori and reappeared only 
after disregarding non-physical degrees of freedom. One of the way of developing the 
standard Hamiltonian description is the symplectic approach by Kijowski and Tulczyiew 
[88], where it is noted that surface integrals are not less important than the volume ones. 
Jezierski and Kijowski developed this approach in GR [89] - [91]. They use an "affine 
formulation", where the connection coefficients F^^ are used rather, than the metric ones. 
The gravitational field is considered inside a closed tube, at a boundary of which some 
conditions are fixed to construct a closed Hamiltonian system. The Hamiltonian describes 
the full energy inside the boundary and has a quasilocal sense [90]. In the linear gravity, the 
requirement of positiveness of the Hamiltonian [89] leads to a "localization" of gravitational 
energy with unique boundary conditions. As an application, gravitational waves on the 
background of the Schwarzschild geometry were studied [91]. 

Based on the symplectic method Nester with co-authors [92] - [98] developed a so- 
called 4-covariant Hamiltonian formulation both for GR and for generaUzed geometrical 
gravitational theories. The Hamiltonian on-shell is a surface integral, which defines a 
quasilocal conserving quantity inside a closed volume. For this approach a displacement 
4-vector constructed from the lapse and shift, and a flat space defined at the boundary of 
the volume are necessary. For the recent development and achievements of this fruitful 
approach see the review paper [97]. 

Returning to the discussion of the previous subsection, we recall again on a necessity to 
operate with local quantities in cosmological and astrophysical applications. Therefore, to 
conclude the subsection let us formulate also the next task: 

• It is important to connect local conserved characteristics with non-local quantities, 
which appear in the theoretical considerations. 

1.3 Goals of the review and plan of the presentation 

Analyzing the problems accented in the previous subsections 1 . 1 and 1.2 a necessity in a 
generalized and universal approach for describing perturbations, both in GR and in gener- 
alized metric gravitational theories, becomes evident. A description, where perturbations 
in a geomettical theory are considered on a curved or flat background, in fact, converts this 
theory into the rank of a field theory, like electrodynamics in a fixed spacetime. A set of all 
the perturbations acquires the sense of the dynamic field configuration. Then, it is desirable 
to represent the perturbed gravitational metric theory in the field-theoretical form (or sim- 
ply, ^eZ<i form) with all the properties of a field theory. We formulate these properties as the 
following requirements: 

(a) The field-theoretical formulation has to be covariant. 



Nonlinear Perturbations and Conservation Laws. 



7 



(b) One can use an arbitrary curved background spacetimes (solutions to GR or another 
metric theory). 

(c) The perturbed system has to be represented as a dynamic field configuration, which 
is associated with Lagrangian and corresponding action. 

(d) The field equations (perturbation equations) have to be derivable from the action 
principle. 

(e) The conserved quantities and conservation laws also have to be derivable using the 
variational principle. 

(f) The field-theoretical formulation has to have gauge freedoms. Gauge transformations 
and their properties have to be connected with the action. 

(g) In order not to have restrictions in the use of orders of perturbations, it is required 
to have an exact formulation for perturbed equations, conservation laws and gauge 
transformations. 

(h) Lastly, it is desirable to have a simple and expUcit form convenient for applications. 

Only such a derivation will permit the required universality and give a full description of 
perturbations. Of course, the field-theoretical formulation has to be equivalent to the geo- 
metrical one, without changing the physical content of the theory. 

Thus, the goal of the paper is to suggest an approach, which gives a possibihty to 
present a perturbed metric (geometrical) gravitational theory in a field-theoretical form, 
which satisfies the above requirements (a) - (h). In last two decades, all the necessary 
elements of such an approach were developed in works by the author together with his 
co-authors, and in other works all of which will be cited later. Therefore, from one point 
of view, the present paper is a review of these works. On the other hand, the paper just 
unites these works into a generalized and universal approach. Our task is to give an outline 
of mathematical development of the approach with necessary mathematical expressions, 
to demonstrate possibilities of the approach and its advantages, and to outline some of its 
applications. Therefore we hope that the present work could be interesting both to the 
experts in gravitational physics, e.g., in conservation laws in GR, and to cosmologists and 
astrophysicists studying the evolution of perturbations on curved backgrounds. This paper 
is not a review of all the numerous perturbation approaches and methods developed during 
the history of GR. Therefore we apologize to the authors whose works are not referred here. 

Let us discuss some of important points. First, a possibility to use an arbitrary curved 
background means that any solution of the initial mettic theory can be considered as a 
background. Thus, the background can be flat, curved vacuum, or even curved including 
background matter, i.e. it can be arbitrary. Second, as was accented, the field and geomet- 
rical formulations of the theory have to be equivalent. This means that a solution of the 
field formulation united together with a background solution have to be transformable into 
a solution of the geometric formulation (initial metric theory). Symbolically this situation 
is explained on the figure 1. Let the slightly sloping curve be related to a background and let 
the oscillating curve mean a solution in the geometrical form, then the difference between 
them symboUzes the solution in the field-theoretical form. Third, usually it is assumed that 
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Figure 1: A symbolical connection between solutions of the same metric theory in the 
geometrical and field-theoretical forms. 

a perturbation of a quantity is less than a quantity itself. We do not impose this restriction 
here. A realization of the requirements (a) - (h) means that the field-theoretical formulation 
can be thought of as an independent exact field theory. Then, of course, the "amplitude" of 
an exact solution of the exact theory can be more than "amplitude" of the background (see, 
e.g., the right side of the figure 1). 

At the earlier stage, Einstein was trying to construct a gravitational theory as a field 
theory in Minkowski space, in the framework of special relativity. However, step by step 
he had concluded that one needs to operate with curved spacetime only, Minkowski back- 
ground space had disappeared from the consideration all together. The construction of 
the field-theoretical formulation in the framework of the geometrical theory is, in a defi- 
nite sense, the revival of the special relativity view. Moreover, as we remarked above, the 
field-theoretical formulation really is an independent field theory and can be constructed by 
independent ways (see below the discussion in subsections 2.1 and 2.2). However, there is 
no contradiction here. The geometrical and field-theoretical formulations are two different 
formulations of the same theory with the same physical content. The background spacetime 
turns out non-observable and has a sense of only an auxiliary structure. 

Up to now the Einstein theory retains its position as the most popular theory of gravity, 
leading in all the applications. Thus, in this paper we pay a significant attention to con- 
structing the field-theoretical formulation for GR which is presented in detail. On the other 
hand, due to the rising precision of experiments and observables in cosmos and due to the 
great interest to the brane models, other metric theories generalizing GR become more and 
more popular and necessary. One of the main properties of the approach presented here is 
its universality. We use this advantage and develop the field theoretical approach applied to 
an arbitrary metric theory. 

The paper is organized as follows. The next section 2 is devoted to the detail descrip- 
tion of the field-theoretical formulation of GR on an arbitrary curved background, which 
have all the properties of a self-dependent field theory. The essential attention is paid to an 
invariance with respect to exact gauge transformations. The last do not effect both coordi- 
nates and background quantities in the field-theoretical formulation and are connected with 
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the general covariance of GR. 

In section 3, we construct conservation laws for perturbations in GR in the framework 
of the field-theoretical approach. Conserved currents and corresponding superpotentials are 
presented. As an important instruments we use the canonical Noether method and the Be- 
linfante symmetrization prescription. The conserved quantities and conservation laws are 
used for examination of asymptotically flat spacetime at spatial infinity, the closed Fried- 
mann model, the Schwarzschild solution and linear perturbations on FRW backgrounds. 

In section 4, we develop the field-theoretical approach to describe perturbations in an 
arbitrary D-dimensional metric theory of gravity on a fixed background. We construct gen- 
eralized conserved currents and corresponding superpotentials with again the essential use 
of the Noether and Belinfante methods. The conserved quantities are tested in the Einstein- 
Gauss-Bonnet gravity for calculating the mass of the Schwarzschild-anti-de Sitter black 
hole. 

1.4 Notations 

Here, we present notations, which will appear more frequently and which are more impor- 
tant. In the text, together with these notations, numerous other notations will be used, they 
will be outlined currently. 

• Greek indexes numerate 4-dimensional spacetime coordinates as well as D- 
dimensional spacetime ones. Usually means a time coordinate, whereas small 
Latin indexes from the middle of alphabet i, j, k, . . . mean 3-dimensional space co- 
ordinates or {D — l)-dimensional hyperspace coordinates; 

• Large Latin indexes A, B, C, ... are used as generalized ones for an arbitrary set of 
tensor densities, for example, = {sf—gg^^, Tap}; 

• A dynamic metric of a metric theory is g^iy {g = det gui/); 

• Bar means that a quantity "Q " is a background one; 

• Thus, g^^ (g = det 5^^) is a background metric. Indexes of all the quantities of a 
perturbed system are raised and lowered by the background metric; 

• Many expressions are presented as densities of the weight -i-l. The reasons are as fol- 
lows. First, all the Noether identities, which are explored intensively, are such densi- 
ties in the initial derivation. Second, conservation laws have to be covariant, however 
partial derivatives are crucial for application of the Gauss theorem. Divergences of 
both vector densities (conserved currents) and antisymmetric tensor densities (super- 
potentials) have this duality property. To accent these expressions we use "hats", as 
more economical notations in this situation. Thus, a quantity "Q^" is such a density, 
it could be a tensor multiplied by -y/^ or -v/^ (for example, Q'^^-^ = \/-^Q°'^'y), 
or could be be independent on these determinants, a situation will be clarified 
from the context; 
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• rjni, is a Minkowskian metric in the Lorentzian coordinates; sometimes is used 
explicitly instead of 1 to stress that a quantity, say 0"^7 = V—vQ'^^j' is a density 
of the weight +1; 

• l^'^ = g^'^ — g^^ = ^—gg'^^ — y/—gg^'^ is a more important form of the metric 
perturbations; 

• Partial derivatives are denoted by (^a), or 9^, 9a,; 

• Da and Da are covariant derivatives with respect to g^i, and with the Chistoffel 
symbols T^^ and T^^, respectively; 

• A^^ = — r'^^ is the tensor actively used in the paper; 

• is an arbitrary displacement vector, whereas A" is a Killing vector of a back- 
ground; 

• The Lie derivative is defined as 

£^Q^ = -eDaQ^ + iDa^^ , 

P 

note the opposite sign to the usual one, Q ^ is defined by the transformation prop- 
erties of Q^; 

• Ctiv = —\ £(gij,v is the tensor actively used in the paper; 

• The Lagrangian derivative is defined as usual: 

uPv, Rfxv: R and R npi,, Rjj,!/, G^i,, Tni,,R are the Riemanman, Ricci, 

Einstein, matter energy-momentum tensors and the curvature scalar for the physical 
and background spacetimes. 

• Usually index "L" means a linearization, for example, G^^ and <I>^^ mean linearized 
pure gravitational and matter parts of the gravitational field equations; 

• The conserved currents are defined in the framework of different approaches as fol- 
lows. In the case of GR: J^^^ is defined with the use of the canonical Noether pro- 
cedure; J^'^^ is defined with the use of the field-theoretical prescription, based on the 
symmetrical energy-momentum tensor; J^'^^ is the canonical current corrected with 
the use of the Belinfante method. The correspondent superpotentials in GR are J'^^y 

• Analogous currents in superpotentials in an arbitrary Z)-dimensional metric theory 
are respectively T'^^y I^^y Jf^) and I^^y l^^^y I^^y 

• (a/3), (ik) and [a/3], [ik] mean symmetrization and antisymmetrization; 



• K — the "Einstein" constant both in GR and in an arbitrary metric theory. 
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2 The exact field-theoretical formulation of GR 

2.1 Development of the field approach 

The study of perturbations in GR, in fact, was begun by Einstein himself. However, as 
a separate field, the history of the field-theoretical approach in GR began in 40's — 50's 
of XX century. Perturbed Einstein equations are rewritten as follows. Define the metric 
perturbations on a flat background in the Lorentzian coordinates as l^^ = g^^ — yf^^tf^; 
l^'^ = {V~v)~^^-^'^ ■ The terms linear in metric perturbations are placed on the left hand side 
of the Einstein equations, whereas all the nonlinear terms are transported to the right hand 
side, and together with a matter energy-momentum tensor are treated as a total (effective) 
energy-momentum tensor t^*^*^. Then Einstein equations are rewritten in the equivalent 
perturbed form as 

G^. = i^ti^f (2.1) 
where, raising the indexes by 77"^, one has the left hand side in the form: 

Gf = ill'^^'^a + Ty'^'^r^.^^ - l-^^'^a - l'"''^,a) ■ (2.2) 

Its divergence identically is equal to zero: duC^ = 0. Then, one obtains directly the 
differential conservation law 

^^^&)=0. (2.3) 

This picture was developed in a form of a Lagrangian based field theory with self- 
interaction in a fixed background spacetime, where is obtained by variation of an 
action with respect to a background metric. Following the introduction in the Deser work 
[99], below, we shall present the main steps in this derivation, the corresponding bibliogra- 
phy can also be found in [99]. Assume that a field theory of gravity in a Minkowski space is 
constructed. By known observable tests (see, e.g., textbook [3]), the most preferable type of 
the gravitational field is the tensor field, say P*^. The linear (approximate) equations have 
to have the form G^^ = and are defined by the quadratic Lagrangian Z^2) • Keeping a 

symmetrical energy-momentum tensor of matter fields as a source of G^^ one obtains 

Gfi,{l) = KT^,ict),n) . (2.4) 

Identically d^Gj^" = 0, therefore Sj/T^ = 0. However, there is a contradiction between 
the conservation law d^Tj^ = and equations of motion for interacting fields cp^. How 
does one avoid this? The right hand side of Eq. (2.4) is to be obtained by variation (con- 
ventionally) both with respect to r]'^'^ and l'^'^. Therefore one needs to make an exchange 
I^A^ j^Mi^l _^ |(/)^, rjf^" + P''! both in the matter Lagrangian and in the matter energy- 
momentum tensor. This just means the universality of gravitational interaction. Next, one 
has to include the gravitational self-interaction. Therefore one adds the symmetrical energy- 
momentum tensor of the gravitational field corresponding to £^2)' ^^e right hand 

side of (2.4) together with T^^((/), r] + I). But the equations that include t^^u^ (l) can be ob- 
tained if a cubic Lagrangian is added, C^^^) + -^(3) • After this one needs to consider the next 
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level, and so on. In the result, one obtains the final variant of the gravitational equations: 



T^,{ct>,l + r,) + J2ti'J'{l) 



n=2 



(2.5) 



It turns out that the equations (2.5) are exactly the Einstein equations, i.e. equations (2.1). 
After the identification ^—rjrj'^'^ + = \/—gg^^ one has only the dynamical metric 
g^^, whereas the background metric r]^^ and the field l^^ completely disappear from the 
consideration. 

Deser himself [99], unlike (2.5), has suggested the field formulation of GR without 
expansions. As dynamical variables he used the two independent tensor fields l^^ and A^^ 
of the 1-st order formalism. After variation of the corresponding action he derives the 
equations in the form: 

G^.(0 = '^[Ca,A) + Cl , (2.6) 



instead of (2.5). After identifications ^/^rj^"" + If^" = ^/^g'^" and A^^ = T^^ the 
equivalence with the Einstein equations is confirmed, but only in the Palatini form, where 
g^^^ and are used as independent variables. 

In the work [100] we have generalized the Deser approach [99]. Instead of the back- 
ground Minkowski space with the Lorenzian coordinates we consider an arbitrary curved 

_ j\ 

background spacetime with a given metric g^j^ and given matter fields $ satisfying the 

background Einstein gravitational and matter equations. We also use the 1-st order formal- 
ism. The gravitational equations get the generahzed form: 

GUI) + ^auil. 4>) = i^ \iUU A) + (2.7) 



where the left hand side is linear in /^'^ and (p^^ and defined later in Eqs. (2.25) - (2.27). The 
term <l>^j^ appears due to <^ . Equivalence with GR in the ordinary derivation is stated after 

identifications: + f^^^ = g^"' , f"^ + A;^^ = T"^^ and ¥^ + (/>^ = In the following 
years we have developed the principles of the work [100] and have used this approach in 
many applications. Our results are presented in the papers [101] - [121], on the basis of 
which, in a more part, the present review was written. 

Elements of the field approach in gravity are also actively developing nowadays in other 
approaches. Thus, in [122], a requirement only of the first derivatives of metrical perturba- 
tions in the total symmetrical energy-momentum tensor has led to a new field formulation 
of GR in Minkowski spacetime, which is different from the formulation in [100]. The new 
total energy-momentum tensor is the source for the non-linear left hand side. On the basis 
of this new field formulation an interesting variant of the gravitational theory with non-zero 
masses of gravitons was developed [123]. A comparison and a connection of the works 
[100, 122, 123] are discussed in [124] in detail. In [125], the work [100] was developed 
to construct the total energies and angular momenta for d + 1-dimensional asymptotically 
anti-de Sitter spacetime. The properties of the field approach [100] appear independently in 
many concrete problems. For example, in [126, 127] a consideration of linear perturbations 
on FRW backgrounds leads to the linear approximation of the exact field formulation of 
GR [100]. In [128], as a development of the field approach, a class of so-called "slightly 
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bimetric" gravitation theories was constructed. In [129, 130], a behaviour of Ught cones 
in Minkowski space and effectively curved spacetimes was examined. Then, based on the 
causality principle, a special criterium was stated. In [131] this criterium was used to show 
that if spatially flat FRW big bang model is considered as a configuration on a flat back- 
ground, then the cosmological singularity is banished to past infinity in Minkowski space. 
The references to earlier works and the theoretical foundation for the field approach can be 
found in the works [132] - [134]. To the best of our knowledge up to date bibhography 
related to the field approach in gravity can be found in [128] - [131]. 

2.2 Various directions in the construction 

There are various possibilities to approach the field-theoretical formulation of GR, which 
are based on different foundations. Here, we shall discuss the well known and important 
ones. The principle used by Deser [99] could be formulated as follows: 

• The source of the linear massless field of spin two (of gravitational field) in 
Minkowski space is to be the total symmetrical (metric) energy-momentum tensor 
of all the dynamical fields, including the gravitational field itself. 

Using this principle Deser has constructed a corresponding Lagrangian and field equations 
and energy-momentum tensor following from it. We have generalized the Deser approach 
on arbitrary curved backgrounds [100]. His principle has been reformulated in a way that 
the linear left hand side of perturbed gravitational equations has to be of the form in Eq. 
(2.7), i.e. together with G^^ one has to include the, Unear in matter perturbations, part 
The analogous principle was suggested for perturbed matter equations. 

The next known method was most clearly presented by Grishchuk [134] and shortly can 
be formulated as: 

• A transformation from gravistatic (Newton law) to gravidynamics, i.e. to a relativistic 
theory of gravitational field (general relativity), equations of which (Einstein equa- 
tions) describe gravitational waves. 

Following this direction, one has to transform the Newton law Ac/) = —AirGp into a special 
relativity description. To satisfy the relativistic requirement a) the mass density p has to be 
generalized to 10 components of the matter stress-energy tensor T^i,; b) the single compo- 
nent (/) should also be replaced by 10 gravitational potentials P'^; c) the Laplace operator 
should be replaced by the d'Alembert operator; d) the gravitational field has to be nonlinear 
and, thus, has to be a source for itself. Following this reformations one obtains the gener- 
alized equations: ^Z^jy'^a = + ^^i/)- They imply that the gauge condition Z'*''^ = 
is already chosen, e) To reconstruct the gauge invariance properties it is necessary to add to 
the left hand side the terms ^ p ~ ^"'ii,v,a ~ ^%/i,a)- ^ result, one obtains Eq. 

(2.1), which is just the Einstein equations. 

The field formulation of GR can be also constructed based on the gauge properties 
(see our work [107]). This direction is analogous to the method in the gauge theories of 
the Yang-Mills type, which are constructed by localizing parameters of a gauge group. 
However, unlike usual, we postulate a non-standard way of localization, namely: 
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• A "localization" of Killing vectors of the background spacetime. 

This assumes the existence of a fixed background spacetime with symmetries presented by a 
Killing vector field A", in which initial dynamic fields (f>'^ are propagated. It is noted that an 
action for the initial fields is invariant, up to a surface term, under the transformation cf)"^ — > 
£x(t)'^. Next, the Killing vector is changed for an arbitrary vector that is localized. Then 
the invariance is destroyed. To restore it the compensating (gauge) field has to be included. 
In doing so the coordinates and the background metric do not change. The requirements to 
have the gauge field as an universal field and to have the simplest sought-for action for the 
free gauge field lead just to the field formulation of GR, developed in [100]. 

As a rule, a fixed background spacetime, in which perturbations are studied, is de- 
termined by the problem under consideration (see subsections 1.1 and 1.2 in Introduction). 
Thus the background could be assumed as a known solution to the Einstein equations. Then 
one has to study the perturbed (with respect to this background) Einstein equations. As al- 
ready was noted above, this picture can be developed as a Lagrangian based field theory, 
where the first step is: 

• The decomposition of dynamical variables of GR into background variables and dy- 
namic perturbations. 

This method is evident in itself and has the explicit connection with the ordinary geomet- 
rical formulation of GR. Namely this method is easily adopted for constructing the field- 
theoretical formulation in the framework of an arbitrary metric theory. In the next subsec- 
tions, basing on the works [100, 103, 104], we present it in detail. However, although the 
construction was developed for both in the 1-st and in the 2-nd order formalisms, here, we 
use the 2-nd order formalism only since it is more convenient and suitable. To the best of 
our knowledge, Bamebey [135] was the first who suggested to use the 2-nd order formaUsm 
for an exact (without expansions) description of perturbations in GR. 

2.3 A dynamical Lagrangian 

Let us consider the Einstein theory and write out its action 




depends on and their derivatives up to a finite order. Then, the gravitational and 
matter equations corresponding to the Lagrangain (2.8) are 



1 SR SC 



0, 



(2.9) 



2k 5gi^'' 5gi^'' 

5C^ 5C^ 



0. 



(2.10) 



The form of Eq. (2.9) corresponds to the form of the Einstein equations 




(2.11) 
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whereas a more customary form is 

{Gui, - nT^y) = G^y - nf^y = , (2.12) 

which is obtained by variation with respect to g^^'^ . 

Next, let us define the metric and matter perturbations with the use of the decomposi- 
tions: 

gi'^ = T + l''\ (2.13) 
$^ = $^ + ^^. (2.14) 

It is assumed that the background quantities g'^" and are known and satisfy the back- 
ground (given) system, which is defined as follows. Its action is 

S= - j d'^xl^ = — - [ d'^xR + - [ d'^xl^ . (2. 15) 

C J ZKC J C J 

The corresponding to the Lagrangian in (2.15) background gravitational and matter equa- 
tions have the form of the barred equations (2.9) and (2.10): 

1 5R Sl^ ^ 
W + 7n^ = 0' (2-16) 



2k 5g^ 5g 

—^ = 0. (2.17) 

Frequently we use a Ricci-flat background with the background equations 

G^. = \. = 0. (2.18) 

Now let us classify the perturbations l^^ and cj)^ as independent dynamic variables, 
which present a field configuration on the background of the system (2.16) and (2.17). To 
describe this dynamical configuration we construct a corresponding Lagrangian called as 
a dynamical one [104]. After substituting the decompositions (2.13) and (2.14) into the 
Lagrangian of the action (2.8) and subtracting zero's and linear in l^^ and (j)^ terms of 
the functional expansion of the Lagrangian C,^ one has 



tdvn ^ t^(^+l^ ^ + ^^E_ Q^f,a ^ _ ^ . (2.19) 

8g (5$ 2k Ik 



As is seen, zero's order term is the background Lagrangian, whereas the linear term is 
proportional to the left hand sides of the background equations (2.16) and (2.17). 

In Eq. (2.19) a vector density is not concreted. However, consider defined as 

k^^r^^-r^l. (2.20) 

with Aj^^ presented as the perturbations of the Cristoffel symbols 

A^, = - = iy"^ i^^gpy D,gp^ - D,g^,) , (2.21) 



16 



A.N. Petrov 



and depending on l'^'^ through the decomposition (2.13). Then a pure gravitational part in 
the Lagrangian (2.19) is presented in the form: 

= -(A^,-A^XW>^ + (r' + ^n(A^.A^<,-A^,A^,) . (2.22) 

It depends only on the first derivatives of the gravitational variables Z^*^. In the case of a 
flat background the Lagrangian (2.22) transfers to the covariant Lagrangian suggested by 
Rosen [136], which has been rediscovered in [137] and [101]. The matter part of (2.19) is 

/ \ /irM ATM 

= £^ (5 + $ + 0) - - - . (2.23) 

2.4 The Einstein equations in the field formulation 

The variation of the action with the Lagrangian (2.19) with respect to l"^ and the contraction 
with -/—ViS^S^ — ^Vni^g"^) give the field equations in the form: 

They coincide with the form (2.7), only now in the 2-nd order formalism. The left hand 
side of Eq. (2.24) is linear in l^^^ and (f)-^. It consists of the pure gravitational part 

= liDfDX. + g^.D^Dj" -D,dJ^ -D,D„i;) , (2.26) 
which is the covariantized expression (2.2), and of the matter part 



It disappears for Ricci-flat backgrounds (2.18), and then Eq. (2.24) acquires the form of 
Eq. (2.1). The right hand side of Eq. (2.24) is the symmetrical energy-momentum tensor 
density 

s:rdyn r) / 1 \ 

f[tot) _ 2— = 2 ( - — t^ -\-t^\=P A-f^ (2 2S) 



The explicit form of the gravitational part is 

V=V {-Sf,S: + hg^.r") (K^^ip - A^^jAf + Dr%J\ ; (2.29) 



+ K [^IP - ^Ppr^da.) + {^Ip - ^%r'9a,) ■ (2.30) 
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The matter part is expressed through the usual matter energy-momentum tensor density T^i, 
of the Einstein theory as 

Taking into account the definitions (2.27), (2.28) and (2.31) in the field equations (2.24) 
one can rewrite them in the form: 

= « (if,, + = i^i^P ; (2.32) 

^^ixu — ^jxv iJ-u '^ixu 29P''^g^ ^pa 2^/^'^^'' i^P'^ 2Sp<^g -^Xr^ T^i, . (2.33) 

The equation (2.32) has the form of Eq. (2.1) even on arbitrary backgrounds. The price is 
that the effective source i^^u^\ including the matter part as Si^^, does not follow from the 
Lagrangian (2.19) directly. However, this matter part could be classified as a perturbation 
of 



= ^ ^— . (2.34) 

Return to the equations (2.24) in the whole. Transfer the energy-momentum tensor density 
to the left hand side and use the definitions (2.25), (2.27) and (2.28) with (2.19): 



■pa 



dg^" SIP'' 



6 / 1 — —M^ 
+ '^^7=nui-7rR + ^ )=0- (2-35) 



Because the second hne contains the operator of the background Einstein equations in (2. 16) 
one can assert that Eq. (2.24) is equivalent to the Einstein equations (2.9). 

By the same way the perturbed matter equations are constructed. They have the form: 

= tA (2.36) 



where 



$^(Z,^) ^ -^l/P- +<^« I, (2.37) 
n = (2.38) 
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2.5 Expansions 

The methods of the exact field-theoretical formulation give a possibility to construct an 
approximate scheme easily and in clear expressions up to an arbitrary order in perturbations. 
Let us show this. Assuming enough smooth functions, expand the Lagrangian jC.^{g+l, 
(p) as 

2! sf^ ^9 5^ 2! ^^A 

+ -T^^r^—^lP''^ + ... + div. (2.39) 
3! 5f sT ^9 

Then the dynamical Lagrangian (2.19) acquires the form: 

First, note that a divergence in (2.40) is not so important because divergences vanish un- 
der the Lagrange derivative. Second, now we can explain why in the linear terms of the 
Lagrangian (2.19) the background equations are not taken into account before variation. 
Indeed, the Lagrangian jC^{g + + in (2.19) contains the same linear terms with 
only the opposite sign not expUcitly (the formula (2.39) shows this). Therefore, in fact, 
these linear terms are compensated, and the real lowest order in the expansion of C'^y^ is a 
quadratic one (2.40). In this relation, recall that in the usual geometrical description of GR 
the attempt to define the energy-momentum tensor through 5C^ /5g^^ leads to zero on the 
solutions of the Einstein equations. In contrast, in (2.28) does not vanish on the field 
equations. The reason is just in the presence in the Lagrangian (2.19) of these linear terms. 

Third, under necessary assumptions the series (2.40) can be interrupted at a correspond- 
ing order. Thus the approximate Lagrangian for the perturbed system can be obtained. Its 
variation gives both approximate field equations and energy-momentum tensor. For ex- 
ample, the quadratic approximation of (2.40) gives a possibihty a) to construct the linear 
equations 



{GUI f) + f)) - i V'-'S' + ^'li = » • 



-$i(Z,</>) ^ ^[[P-^ + 0S^| =0, (2.42) 



b) to construct the quadratic energy-momentum tensor: 



(2.43) 
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The cubic approximation of (2.40) gives a possibility a) to construct the field equations 
including quadratic terms (which are related to the energy-momentum tensor), and b) to 
construct the energy-momentum tensor, including quadratic and cubic parts, etc. 

Fourth, the expansions, like (2.39) - (2.43), are used in quantum field theories [4] and 
called as functional expansions. As is seen, in the framework of the classic theory the func- 
tional expansions (2.40) - (2.42) give, in fact, the algorithm for constructing the approximate 
systems, thus an each order can be obtained automatically. 

2.6 Gauge transformations and their properties 

The important property of the field formulation is the gauge invariance. Usually gauge 
transformations in GR and other metric theories are connected with mapping a spacetime 
onto itself that is connected with differentiable coordinate transformations 

x"^ = /"(x^) . (2.44) 
These transformations can be connected with the smooth vector field 

To map the spacetime onto itself one has to follow the standard way [138] - [140]. After the 
coordinate transformations (2.44) (or (2.45)) the metric density, for example, is transformed 
as g^'^ix) ^ g''^'^{x'). Then return to the points with quantities x within a new frame {x'}. 
After that one has to compare geometrical objects of the initial spacetime and of the mapped 
spacetime in the points with quantities x. The comparison can be carried out both without 
in the terms of (2.44) and with included in (2.45): 

oo 

g"^^{x) = gi^^{x) + 6fg'^^=gt^^ix) + J2-£lt^{x), (2.46) 

fc=i 

oo 

^'^(x) = $^(x) + 5/$^(a;) = $^(a;) + ^-4$^(x). (2.47) 

k=i 

The next property is very useful. Assume that geometrical objects ^'^ are differentiable 
functions of other geometrical objects ijj^ and their derivatives, but are not expUcit functions 
of coordinates. Then it is clear that a simple substitution gives '^^{tp'^) = and one 
has 

oo 1 

^^{tP'^ix)) = ^^i^^ix)) + 6f^^ = ^^{tP^{x)) + ri"^?*"*^ • (2-48) 

k=i ^■ 

Now let us define the gauge transformations for the dynamical variables in the frame- 
work of the field formulation of GR: 

oo 1 

f/M^ ^ [M^ + (^1^^ + lUi^^j ^ fM^ + ^ i. £| {^gf^^ + ^ (2.49) 

k=l 

oo -■ 

0'^ = <^^ + ^^($^ + <^^)=<^^ + ^i.4($^ + <^A^. (2.50) 

fe=l 
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The assumptions above, in fact, state that the operators S / and Ylk^i h equivalent. 
However, the operator 5f could present a wider class of transformations, which cannot be 
expressed through infinite series. We conserve a possibility to consider such transforma- 
tions, however more frequently we will use the sums, keeping in mind that they can be 
changed hy 6f also. 

Now, let us substitute (2.49) and (2.50) into the dynamical Lagrangian (2.19). One finds 
that this substitution into jC.^{g + l,^ + (j)) just permits to use the property (2.48). Then 
finally one obtains 

oo 1 



k=l 



5C^ 



■fiv 



2k 



(2.51) 



Because is the scalar density of the weight +1 all the terms under the sum in (2.51) are 
divergences. Thus C'^^^ is gauge invariant up to a divergence if the background equations 
(2.16) and (2.17) hold. 

Considering the gauge invariance properties of the field equations we use their form 
(2.35). The substitution of (2.49) and (2.50) with the use of the property (2.48) gives 



+ 



dg' 



f—pT OO 



Sn 



Sn 



nt 



(tot) 

6-K 



2.^ 

Sg 



(2.52) 



Thus the field equations are gauge invariant on solutions of themselves and with using the 
background equations (2.16) and (2.17). Analogous transformations could be presented for 
the matter equations (2.36). 

Concerning the energy-momentum tensor, it is not gauge invariant. Even on the dy- 
namical equations, as it follows from (2.52), under the transformations (2.49) and (2.50) 
one has 

i^^iT^ = >^itf + G]:,{1' -l) + $^,(Z' -l;<t>'-ct>). (2.53) 

The mathematical reason is that the background equations in the Lagrangian (2.51) cannot 
be taken into account before variation. In the case of the Ricci-flat backgrounds (2.18) 
one has <l^j^ = 0, therefore the energy-momentum is gauge invariant up to G^^^ — 
covariant divergence (see (2.26)). Let us turn also to the equivalent form (2.32) of the 
field equations with the effective source i]j^v^\ For the operator of the field equations 
^]lv ~ i^i^^^v^^ the form of the transformations (2.52) can be used without changing. Then 
on the field equations one has 



(2.54) 



i.e. for all the kinds of backgrounds t}f/'^^ is gauge invariant up to a covariant divergence. 
It is not surprising that both the energy-momentum tensors are not gauge invariant. It is a 
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(b) 



Figure 2: The perturbations (a) and (b) connected by the gauge transformations. 

waited result. Indeed, this reflects the fact that energy and other conserved quantities in GR 
are not localized. Moreover, the formulae (2.53) and (2.54) are very useful because they 
give a quantitative and constructive description of the non-localization. See the discussion 
in subsection 1.2. 

The transformations (2.49) and (2.50) are directly connected with a way of mapping a 
perturbed spacetime onto a given background spacetime (a perturbed solution onto a given 
solution). In the other words, they are connected with a definition of perturbations with 
respect to a given background. Let us consider a solution in the geometrical form g^^{x). 
Next let us map a spacetime onto itself following the prescription at the beginning this 
subsection. Then we obtain g'^^^{x). After that we decompose both of the solutions into 
dynamic and background parts: 

r^ix) = T"{x)+i'"'{x), (2.55) 
(g'^^nix) = T{x) + i'^^x). (2.56) 

The main property of (2.55) and (2.56) is that the background metric is the same. Then 
it turns out that l^^{x) and l'^^{x) are connected by the transformations in (2.49). This 
situation is interpreted as follows. The same background in (2.55) and (2.56) is chosen 
by different ways that symbolically is illustrated on the figure 2. The curves in both the 
cases (a) and (b) symbolize the same solution of GR in the geometrical form, whereas the 
straight lines present a background, say, a Minkowski space. The perturbations in the cases 
(a) and (b) are different, but they are solutions to the equations of the field formulation 
connected by the gauge transformations, and, thus, they are the same solution in different 
forms. In spite of that the gauge transformations in the field formulation have the evident 
geometrical origin, they can be interpreted as inner gauge transformations, like in standard 
gauge theories. Indeed, they act only onto the dynamical variables (perturbations), whereas 
the backgrounds variables and coordinates do not change. 
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As discussed in Introduction, in many of applications it is important to consider equa- 
tions and gauge transformations in linear and quadratic approximations. Assume that per- 



turbations and their derivatives are small: < f, cj)^ < « dj^"" ■ ■ ■ and 



Assuming that the background equations (2.16) give a connection 



g h{K)^ with a coefficient h{K) depending on the Einstein constant one has to set 
l^'^ K, h{K,)(f>'^, etc. Now, rewrite the equations (2.24) up to the second order in perturba- 
tions: 



Assuming iterations the perturbations can be expanded as = If^^ + + 
(f)'^ = (pf + (f)2 + ■■■ ■ Then the linear equations acquire the form: 

GUh) + ^^p{ii, 0i) = O, 
whereas the quadratic equations have the form: 



(2.57) 
. , and 

(2.58) 
(2.59) 



Besides, assuming = + ~^ ■■■ with 9aCi 



jfMU 



h{n)<Pi and 



?2 ~ ^0^2 ~ ■ • • ~ ^2*^ ~ h{K)(j)2, one has a linear version of (2.49) in the simple form: 



Under these transformations the equations (2.58) are transformed as 



(2.60) 
(2.61) 



Gi 



+ 



R 



pa 1^ per 29pcr'^ 



(2.62) 



and are, thus, gauge invariant on the background equations. In the simple case of the Ricci- 
flat background (2.18) the linear transformations have only the form (2.60), without (2.61). 
Then the formula (2.62) transfers to the formula = G^jj, which expresses the gauge 
invariance of the linear field of spin 2. It is the well known gauge transformations in the 
linear gravity [1,3]. 

The quadratic order of the gauge transformations has a form; 
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(2.63) 



Substitution of Eqs. (2.60) and (2.63) into (2.59) give 



+ 



+ 



dg 



■pa 



(^£^2 + ^£li^ Rpa-n(^ 



dgp' 



17% 



Rpa l^{Tpa 29pa'^] 



+ 



(2.64) 
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Thus, equations (2.59) are gauge invariant on the background equations (2.16) and on the 
Unear equations (2.58). Of course, the procedure can be continued in the next orders. 

In this subsection, we were based on the exact theory of gauge transformations devel- 
oped in our works [100, 104]. Together with this, for the presentation we have used some 
of details given in [101, 112, 113]. An arising interest to cosmological perturbations stim- 
ulates their more detailed study including a second order approximation [14]. Thus it turns 
out necessary to examine the gauge transformations up to a second order also. Such stud- 
ies (independently on ours, but repeating them in main properties) were carried out, for 
examples, in [141, 142]. 



2.7 Differential conservation laws on special backgrounds 



The energy-momentum tensor is the one of important objects of a field theory in Minkowski 
space. Its differential conservation together with symmetries of the Minkowski space permit 
to construct integral conserved quantities (see subsection 3.2). The field formulation of GR 
has also the conserved energy-momentum in Minkowski space with the same properties 
(2.3). In this short subsection we demonstrate that the conservation law, like (2.3), also 
has a place on curved backgrounds, which are important in many applications. Although 
conservation laws and conserved quantities are constructed and studied in the next section 
3 in detail, by the logic of the presentation we include this subsection here. 

Firstly, consider Ricci-flat (including flat) backgrounds (2.18), which have an indepen- 
dent meaning. This means that = and = 0. Then the Lagrangian (2.19) is 
simplified to 

1 



Idyn ^ £9 +t^ = £.3 + t 



M 



2k 2k 
The field equations (2.24) transform into the form 



G 



f3 A- f 



(2.65) 



(2.66) 



For (2.18) one has identically Dj,G^^ = and, thus, a divergence of Eq. (2.66) leads to the 
differential conservation law: 

D^tjf"*)'^ = . (2.67) 

Now consider backgrounds presented by Einstein spaces in Petrov's definition [143], 
then the background equations are 



(2.68) 



where A is a constant. Ricci-flat and AdS backgrounds are particular cases of Einstein's 
spaces. The Lagrangian of the background system has the form: 



= R + C^ 

2k 



1 

'2k 



R-2Ay 



-9 



(2.69) 



Here, the constant A is rather interpreted as 'degenerated' matter. Then, the dynamical 
Lagrangian (2.19) transforms into 



Idyn ^ Ig _^ ^r, 

2k 



5g 



(2.70) 
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and leads to the field equations (2.24) in the simple form: 

G^, + A/>^ = K (tl, + ) ^ 4*;*) . (2.71) 
Taking into account the background equations (2.68) one has identically 

D,(g]:^ + M^^=Q. (2.72) 

Thus, the differentiation of Eq. (2.71) gives the same conservation law (2.67). In heuristic 
form the differential conservation law on AdS and de Sitter (dS) backgrounds was used in 
[68]; in the Lagrangian description it was shortly noted in [100]; and, in the paper [144], it 
was studied in more detail. 

2.8 Different definitions of metric perturbations 

In GR, the metric perturbations can be defined by different decompositions: 

r = T + i'"', 

... = ... + ..., (2.73) 
not only by (2.13). Denoting components of metrical densities in the united way 

5" = 9^.u, V^g"", ^99^.. {-9)9^\ ■ ■ •} (2-74) 

we rewrite the action of GR as 

S = - \ d^xt^^"^ = ^ / d^xRig'') + - / d^xL^(^^, g") . (2.75) 

c J 2kc J c J 

After its variation the Einstein equations take the generahzed form: 

1 SR 5C^ ^ ^ ^ 
1^ + ^^ — = 0. (2.76) 

2k 5g"- Sg"- 

The background action and equations have the corresponding to (2.75) and (2.76) barred 
form. After that we present the united form for decompositions (2.73) 

g''=r + h\ (2.77) 

and, following the rules of (2.19), construct the generalized dynamical Lagrangian: 



Nonlinear Perturbations and Conservation Laws. 



25 



The total symmetrical energy-momentum tensor density is defined as usual: 

4*°* ^ 2^ . (2.79) 

After substituting the expression (2.78) into this definition (identity) and taking into account 
Eq. (2.76) and the barred Eq. (2.76) we obtain the Einstein equations in the form (2.24): 

G'm^"^ + ^M^"^ = • (2-80) 

Here, the hnear left hand side is defined by the same operators (2.25) - (2.27), only now 
with independent variables 

rs—lJ.V 

fM- ^ (2.81) 

However, a choice of two different arbitrary decompositions as = 5? + and 
92 =52 + ^2' gives the difference 

'(a2) \a\) ~ P(a)\2 ' l,Z.OZ; 

which is not less than quadratic in perturbations. Because differences enter the linear ex- 
pressions of equations (2.80) the energy-momentum tensor densities fj^"* and 
have the same differences too. For the case of flat backgrounds this ambiguity was noted 
by Boulware and Deser [145]. Later we [104] have examined it for arbitrary curved back- 
grounds and arbitrary metric theories. However, only in our works [1 16, 118] this ambiguity 
has been resolved, and we present this solution in subsection 3.5. 



3 Conservation laws in GR 

3.1 Classical pseudotensors and superpotentials 

During many decades after constructing GR pseudotensors and superpotentials were main 
objects in constructing conservation laws and conserved quantities. In the framework of the 
field-theoretical approach these notions and quantities, in a definite sense, are developed 
and generalized. Therefore in this subsection we give a short review describing classical 
pseudotensors and superpotentials, only which are necessary in our own presentation. On 
their examples we outUne the general properties of such objects, their problems and some 
of modem appUcations. 

Let us present the general way for constructing pseudotensors and superpotentials. Us- 
ing the metric and its derivations consttuct an arbitrary quantity lA^'^, satisfying the condi- 
tion 

a^aWr^O. (3.1) 

Next, define the complex 

= daUr - -G^ , (3-2) 
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which usually is called as an energy-momentum pseudotensor of gravitational field. Using 
the Einstein equations one obtains 

+ 0^ = daUr (3-3) 

where plays the role of a superpotential. Thus Eq. (3.3) is another form of the Einstein 
equations. Due to (3.1) one has a differential conservation law 

d^{fl^ + 0i:)=O. (3.4) 

Concerning the physical sense, it is the 4-dimensional continuity equation, which is the 
differential conservation law derived directly from the field equations. As is seen, the above 
construction has problems. First, it is an ambiguity in a definition of 6^ and iY^", second, 
these expressions are not covariant in general. 

The presented picture, as a generaUzation, has been formulated after prolonged history 
of constructing pseudotensors and superpotentials in GR (see a review [146]). Below we re- 
call constructing some of well known expressions. In the work [55], following the standard 
rules Einstein had suggested his famous pseudotensor Firstly, a non-covariant Lagrangian 



^cut _ 

~ 2k 



— LflA*!^ (r'^ TP - V T"' ) (3 5) 



had been suggested. It differs from the covariant Hilbert Lagrangian 

t" = -—R (3.6) 

by a divergence and leads to the same field equations (2.12). Next, the corresponding to 
(3.5) canonical complex had been constructed: 

pi r cut 

which is just the Einstein pseudotensor Both the Lagrangian (3.5) and the pseudotensor 
(3.7) depend on the metric g^j^i, and only its first derivatives. Combining (3.7) and the field 
equations Einstein had obtained the conservation law 

d^(fl^ + i^^') =0, (3.8) 

which is a variant of (3.4). 

Later Tolman [147] had found the quantity 7^^" with a property (3.1), and for which, 
as a particular case of (3.3), one has Tif + = da'I'Jf''^- From here the conservation 
law (3.8) follows directly. However, T^"' is not antisymmetrical in a and /3, whereas the 
use of antisymmetrical superpotentials = —U^'^ is more reasonable because it makes 
expUcit the identity (3.1) and presents less difficulties under covariantization, etc. Assuming 
an antisymmetry Freud [148] had found out the superpotential 



2k V-5 



(-5)(rt"^-5^V^)] , (3.9) 
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such that 

fj}+i^i' = daPr, (3.10) 

and connected with the Tolman one by = + S/jd^^"^; 5„04>{^°^ = 0. 

Many authors (see, e.g., more known publications [149] - [151] and the review [146]) 
considered the problem of an ambiguity in definition of conserved quantities described in 
(3.1) - (3.4). It has turned out that some of pseudotensors and superpotentials are directly 
connected with a Lagrangian through the Noether procedure. This restricts some ambigui- 
ties in their definitions. Naturally, the Eienstein pseudotensor is defined by the Lagrangian 
(3.5) which, being generally non-covariant, is co variant with respect to linear coordinate 
transformations. Then for translations X^^-^ = 5^ one can write out the Noether identity: 

£^t^' + d^(x^^^^t^')=0. (3.11) 

Direct calculations lead to + i^^^ = with the definition (3.7) that after us- 

ing Eq. (2.12) gives the conservation law (3.8). The general conclusion is formulated as 
follows. 

• In the sense of the Noether procedure the pseudotensor (3.7) is uniquely defined by 
the Lagrangian (3.5); other Lagrangians C = Z!'^"* + div lead to other pseudotensors. 

Thus the identity (3.11) with the Hilbert Lagrangian (3.6) leads to 

instead of (3.10), where the M0ller [151] superpotential 

xr = ^^99^^ 9"^' [d0.p - dp9.p) (3.13) 

is used. The components M.^[x°') of the pseudotensor in Eq. (3.12) are transformed as a 
4-dimensional vector density under transformations x'^ = f^{x^),x'^ = x^ + g{x^) [151]. 
In contrary, the corresponding components of the Einstein pseudotensor (3.7) do not. This 
property of 7V4[^(x"), in apart, improves the general non-covariant picture. 

The other problem is in a definition of angular momentum of a system. Both Einsten's 
and Mceller's pseudotensors themselves, as canonical Noether complexes, cannot define it. 
One needs to modify these expressions or to construct new ones (see the review [146]). 
One of methods to construct angular momentum, using unique energy-momentum com- 
plex, is to construct symmetrical pseudotensors. Thus, Landau and Lifshitz have suggested 
a such expression [1]. Their famous symmetrical pseudotensor t^^, like i^*^, has only the 
first derivatives, and is expressed through Einstein's (3.7) and the Freud's (3.9) expressions: 
*LL ~ g^pT^'^^ + g'^Hp'^. However, as is seen, it has an anomalous weight +2, un- 
like expressions in (3.10) and (3.12) with the weight +1. Goldberg [150] generalized the 
Landau-Lifshitz approach, suggesting a symmetrical pseudotensor of an arbitrary weight. 

Concerning symmetrical expressions, return also to the the equation (2.1) and the ex- 
pression (2.2). The last can be rewritten as 

Gf = Kdpf''^'''^^ (3.14) 
ft^i'^P] = ig^ (^[^'[m^"]/^ _ . (3.15) 
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One can find that is the well known superpotential by Papapetrou [152]; in the Unear 

approximation it is Weinberg's superpotential [153]; and, up to factor ^/^, it coincides 
also with the linear version of the Landau-Lifshitz superpotential [1]. Then the equation 
(2.1) can be rewritten as 

This means that the conservation law (2.3) can be interpreted by the same way as the con- 
servation law (3.4), i.e. in the framework of the superpotential approach. 

It was a problem to relate modified superpotentials to Lagrangians, like it was done for 
Freud's (3.9) and M0ller's (3.13) ones. Only recently with the use of a covariant Hamil- 
tonian formalism [94] (see subsection 1.2) Chang, Nester and Chen [95] closed this gap. 
They describe GR by Hamiltonians with different surface terms corresponding to different 
boundary conditions under variation. For each of known (covariantized) superpotentials 
one can find its own Hamiltonian. Thus, the Dirichlet conditions correspond to the Freud 
superpotential (3.9), Neuman's ones — to M0ller's superpotential (3.13). All the other 
superpotentials require more complicated boundary conditions. Recently Babak and Gr- 
ishchuk [122] suggested a special Lagrangian including additional terms with Lagrangian 
multipliers, which leads to a covariantized Landau-Lifshitz pseudotensor of normal weight 
+L 

Of course, a non-covariance of pseudotensors and superpotentials lead to the evident 
problems described in all the textbooks. In this relation, return to the Hilbert Lagrangian 
(3.6) recalling that it is generally covariant. Thus, for constructing conservation laws one 
can use arbitrary diplacement vectors not just the linear translations Aj*^^, like in (3.11). 
Hence there has to exist a generally covariant energy-momentum complex and superpoten- 
tial. Such quantities were found by Komar [154], his superpotential is 

= (L)/^^" - D"^/*) = y^D^i^^"^ = -D^'^i'^^ . (3.17) 

With ^'^ = const it goes to the M0ller superpotential (3.13). However, in spite of the 
advantage with covariance, there are problems. The superpotential (3.17) does not make a 
correct ratio of angular momentum to mass in the Kerr solution (see, e.g., [137, 140]). 

The mentioned above problems could be neglected if resonable assumptions are made. 
Then many classical pseudotensors and superpotentials could be useful for a description of 
differerent physical systems. In this sense it is useful to return to Einstein's arguments. First, 
in spite of that the local conservation laws (3.8) are not covariant they have the same form 
in every coordinate system. Being the continuity equations they give a balance between a 
loss and a gain of a density of particularly defined physical quantities, and hence they state 
the equilibrium of a system. The second remarkable Einstein's example [59] is related to 
masses connected by a rigid rod. Tensions in the road (which exist because is present in 
(3.8)) can be interpreted only as a compensation to tensions of gravitational field included 
in tfj,^'^. Third, Einstein had suggested also a simple method for the "localization" [58], 
when an isolated system is considered as placed into "Galileian space", and "GaUleian 
coordinates" are used for description. Finally, he had used his energy-momentum complex 
to describe weak gravitational waves as metric perturbations with respect to a Mikowski 
metric [56] - [58]. 
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At reasonable assumptions, pseudotensors and superpotentials are left useful nowa- 
days too. Thus, usually quantities of global energy-momentum and angular momentum 
for asymptotically flat spacetimes are calculated. For example, Bergmann-Tomson's and 
Landau-Lifshitz's complexes give the same expressions for global angular momentum in 
an isolated system [155]. In [156, 157], with the use of the Einstein, Tolman and Landau- 
Lifshitz complexes reasonable energy and momentum densities of cylindrical gravitational 
waves were obtained. Recently, in the work [158] with a careful set of assumptions it was 
shown that a tidal work in GR calculated for different non-covariant energy-momentum 
complexes is the same and unambiguous. Authors frequently use pseudotensors and super- 
potentials to describe many popular exact solutions of GR (see [159] - [166] and references 
there in). Particulary, it was shown that in many cases different pseudotensors give the same 
energy distribution for the same solution. 

However, in spite of some successes in applying pseudotensors and superpotentials, 
it is clear that more universal and unambiguous quantities are more preferable. The next 
subsection is devoted to possible ways in improving properties of the classical complexes. 



3.2 Requirements for constructing conservation laws on curved backgrounds 

From the beginning we recall how conservation laws are constructed in a field theory in a 
Minkowski space. These principles are useful both to understand how classical pseudoten- 
sors and superpotentials could be modified to avoid their problems and to better describe 
and explain conservation laws in the field-theoretical formulation. 

Consider an arbitrary theory of fields in the Minkowski space in curved coordinates 
with the Lagrangian 



. (3.18) 



ifc). = -j^^D^V - ■ (3.20) 



For the sake of simplicity assume that it contains the first derivatives only. The symmetrical 
(metric) energy-momentum tensor density is defined as usual: 

Hs)a/3 = 2^=^ . (3.19) 

The density of the canonical energy-momentum tensor has the form: 

Both of them are differentially conserved 

DJl^, = 0, (3.21) 
= 0, (3.22) 

on the field equations. To construct integral conserved quantities it is necessary to use 
Killing vectors of the Minkowski space. Contracting i"^^^ with an each of the 10 Killing 
vectors one obtains the currents 

A't)W = ^T^).^^ (3.23) 
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which are conserved differentially 

Djf,)(A)^aJ(",)(A) = 0. (3.24) 

The other situation with the canonical energy-momentum f^f.^^- If one tries to conttact 
it with the 4-rotation Kilting vectors, then the correspondent currents, like (3.23), are not 
conserved. The reason is that i^"^ is not symmetrical in general. A one of the ways to 
construct differentially conserved currents is 

(A) = %).>^'' + ^l^^.DpX^ (3.25) 



where 



(3.26) 



is a spin density. Another way follows the Belinfante symmetrization [167], who has intto- 
duced the antisymmettic in a and (5 expression 

■^(V) -^(i>) +^{i>) -^(v) ^^■^'> 
and has defined the new (symmetrized) energy-momentum tensor density as follows 

iW = ifc). + Mty (3.28) 

It is also conserved, D^i^^^^ = 0, on the equations of motion. Besides, it is symmetrical 
on the equations of motion. Then the symmetrized current 

= (3-29) 

is also conserved for all the Killing vectors, i.e. one can use the unique object t'(^B)(T without 
additional quantities, like the spin tensor in (3.25). 
For the Minkowski background it is easy to show 

%)a = + = iiB)a (3-30) 

on the field equations. Thus the metric energy-momentum (3.19) is equivalent to the Be- 
finfante symmetrized quantity (3.28) (in this relation see also [168]), i.e. = 7^^. Note 

also that I^^^ differs from I^^^ by a divergence. One can see that the Belinfante proce- 
dure suppresses the spin term in the current (3.25). Classical electrodynamics has just the 
Lagrangian of the type (3.18), and is a good illustration of the above. 

The first example where the BeUnfante method was used in GR is the paper by Pa- 
papertrou [152]. He has symmetrized the Einstein pseudotensor and obtained his remark- 
able superpotential. By the original definition [167] it is necessary to use a background 
metric. Thus, Papapetrou used an auxifiary Minkowski metric. Later Berezin [169], in 
the framework of the field approach in GR on a flat background, has shown that effective 
energy-momentum tensor can be constructed by the Belinfante method. Applications of 
the Belinfante method in GR to the pseudotensors without using an auxiliary background 
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Figure 3: 4-dimensional volume V4 restricted by a truncated cylinder. 

metric [168, 170] lead uniquely to the Einstein tensor. It is important and interesting result. 
However, this cannot be useful for description of energy characteristics, e.g., in vacuum 
case. Recently Borokhov [171] generalized the Belinfante method for an arbitrary field 
theory on arbitrary backgrounds, with arbitrary Killing vectors and with non-minimal cou- 
pling. With simplification to minimal coupling his model goes to the standard description 
[60, 168, 170]. 

Denoting Ij^,^^ = 1 ^(^) } consider 4-dimensional volume V4 in the Mrkowski space, 
the boundary of which consists of timelike "wall" S (cylinder) and two spacelike sections: 
^0 •= ^0 = const and Si := ti = const (see figure 3). Because the conservation laws, like 
(3.24), are presented by the scalar densities, they can be integrated through the 4-volume: 
Jy^ dfilfli;^ (A)(i^x = . The generalized Gauss theorem gives 

/ /(%(A)d3x-/ if^^{X)d'x+<f I^^^{X)dS, = (3.31) 
where dSfj, is the element of integration on S. If in Eq. (3.31) 

i i^^^{X)dSf, = , (3.32) 

'J S 

then the quantity 

P(^)(A) = Jj^^)iX)d'x (3.33) 

is conserved on S restricted by dT, — intersection of S with S. If the equality (3.32) does 
not hold, then Eq. (3.31) describes changing the quantity (3.33), i.e. its flux trough 9S. 
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Due to the difference between 1^^^ and /^'^^ the quantities 7^(5) (A) and 7^(c) (A) differ one 
from another by surface integrals. However, as a rule, boundary conditions suppress these 
surface terms, and then 'P(s)(A) = V(^^^{\). 

Now let us return to pseudotenors and superpotentials. To avoid some problems they 
have to be covariantized. Usually a covariantization is carried out by including an auxil- 
iary metric of a Minkowski space. After that partial derivatives are related to Lorentzian 
coordinates of this flat spacetime. In arbitrary coordinates, partial derivatives naturally 
transform into covariant ones. Thus, one can keep in mind that in usual expressions, like 
(3.3) in general form, already there exist the Minkowski metric r/^^/ and its determinant 
r] = det7?^jy = —1. For constructing integral conservation laws based on covariantized 
pseudotensors it is also natural to use the Killing vectors of the Minkowski space. 

Let us present the program of the covariantization in detail. Consider the generalized 
conservation law (3.4) and rewrite it in the equivalent form: 



ri{{0^p+fp/V^)\1,)]=0. (3.34) 



Here, the left hand side is a scalar density. The translation Killing vector A^^^ of the 
background Minkowski space is expressed in the Lorentzian coordinates, therefore one 
has Aj'^^ = 6^ (the lower index is not coordinate one). Thus, in (3.34) the differentially 
conserved vector density (current) 

i'^(A) ^ V=^((^M + fM)/^)Af^), (3.35) 
Df,J^'{X) = df,J^'{X) = (3.36) 

is presented. Here, one has to accent that if 6'^'^ is symmetrical, then the conservation law 
(3.36) could be extended for 4-rotation KilUng vectors. However, if 6^'^ is not symmetrical, 
then the simple and direct covariantization is restricted to using only translation KilUng 
vectors. 

Next, as is seen, the programm of constructing integral quantities (3.31) - (3.33) can be 
repeated for the covariantized current (3.35). The superpotential in Eq. (3.3) can also be 
rewritten in the covariantized form: 

Jf^-{X) = ^ (uj^yV^) Af,) . (3.37) 
Thus the expression (3.3) can be rewritten in a fully covariant form: 

J^(A) = J^"(A) = :d„ J^«(A) (3.38) 

that has a sense of the conservation law because (A) is an antisymmetric tensor density. 
Due to antisymmetry of the superpotential (3.37), the 4-momentum presented in the volume 
form, Uke (3.33), is expressed as a surface integral 

r{X) = i J^^{X)dsk (3.39) 

JdT, 



where dsk is the element of integration on (JS. 
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Up to now in constructing conservation laws we considered Killing vectors of a back- 
ground only. However, other displacement vectors can be used also. Let us give two ex- 
amples. First, for study of perturbations on FRW backgrounds a so-called notion of an 
integral constraint, which connects a volume integral of only matter perturbations with a 
surface integral of only metric perturbations, frequently is very important [172]. With the 
use of such constraints, e.g., measurable effects of the cosmic background radiation were 
analyzed [173]. In the definition of integral constraints integral constraint vectors, not 
necessarily Killing vectors, play a crucial role [172]. Second, in [174], a new conserved 
energy-momentum pseudotensor was found and used in an effort to integrate Einstein's 
equations with scalar perturbations and topological defects on FRW backgrounds for a sign 
of spatial curvature /c = 0. In [175] it was realized that these conservation laws might be 
associated with the conformal Killing vector of time translation, and can be also generalized 
for A; = ±1. 

Taking into account the above we formulate the generalized requirements for construct- 
ing conservation laws for perturbations on curved backgrounds as follows. 

• Expressions have to be covariant on a chosen background and valid for arbitrary 
curved backgrounds. 

• One has to construct differentially conserved currents (vector densities) such 
that d^J^^{Cj = 0, and which are defined by the canonical Noether procedure apphed 
to a Lagrangian of a perturbed system. 

• The currents have to be expressed through corresponding superpotentials J^'^ (an- 
tisymmetric tensor densities), such that dfj,i,J^^{^) = 0, by the way J^{C) = 



• There has to be a possibility to use arbitrary displacement vectors not only Killing 
vectors of the background. 

• Applications of suggested conserved quantities and conservation laws have to satisfy 
the known tests (see discussion in subsection 1.2). 

3.3 The Katz, Bicak and Lynden-Bell conservation laws 

Katz, Bicak and Lynden-Bell [176] (later we call as KBL) have satisfied the requirements 
at the end of previous subsection. They describe a perturbed system by the Lagrangian: 



see notations in (2.8) and (2.15). In fact they use bimetric formalism, where g^^ and 5^,^ 
are thought as independent metric coefficients. The gravitational part of (3.40) is 



(3.40) 




(3.41) 
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It is also very important to note that is defined in (2.20). 

For the Lagrangian (3.41), as for a scalar density, one has the identity: 



(3.43) 



It is a generalization of Eq. (3.11) for arbitrary displacement vectors Consequent iden- 
tical transformations give: 



V = 0, 



(3.44) 



SjCg 



•pa 



^9^ 



pa 



pa 



u 

(c) 



(3.45) 



Note that here 



SCg 

1 pa 



a 1 —^l 



(3.46) 



Then substituting the Einstein equations (2.12) and their barred version into Eq. (3.44) one 
obtains the differential conservation law 



0, 



The first term in (3.48), the generalized total energy-momentum tensor density, is 



(3.47) 
(3.48) 



0(c). = + (Tlf - tT) + Y/'Rpa^i , 



(3.49) 



with lf° 



9' 



pa 



g^'^ . The first term in Eq. (3.49) is the canonical energy-momentum 



tensor density of the gravitational field: 

- r (a^,A^, - A^,A^,) 5>i + (a-,A^, - A^,A^,) , (3.50) 

the second term in Eq. (3.49) is a perturbation of the matter energy-momentum, the third 
term describes an interaction with the background. The second term in (3.48) is the spin 
tensor density: 



2Kai'i"' = - 2k- 



dCc 



(3.51) 
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If one simplifies to a Minkowski background then it is the quantity presented by Papapetrou 
[152] to construct angular momentum with the use of the Einstein pseudotensor. The last 
term in (3.48) has the form: 

^fe)(0 = ^ [i^'dxCp + /"'^ - 2D,C!i)] , (3.52) 

with 2(po- = —£(gpa = 2D(^p^„y It, thus, disappears if .^^ = is a KiUing vector 
of the background spacetime. Besides, in this case the part of the second term in (3.48): 
(T^'^'^L'(pAct) disappears also. 

Because Eq. (3.44) is the identity then the quantity must be presented through a 
superpotential «^(^) (0' ^^^s indeed 

m) = dJ^^^^iO, (3.53) 
Jf^Ji) ^ + (3.54) 

This superpotential has been found earlier independently by Katz for flat backgrounds [137] 
and by Chrusciel for Ricci-flat backgrounds [177]. However, as it turns out, it has the same 
form (3.54) for generally curved backgrounds. Note that the expression (3.54) modifies 
the Komar superpotential (3.17). With the use of the Einstein equations, the identity (3.53) 
transforms into 

Jf^)(0=a,Jf^^(e), (3.55) 

that is another form of the conservation law (3.47). 

Now, let us compare the KBL expressions with the Einstein prescription. For the sim- 
plification (jfiij rjfiiy the KBL Lagrangian jCq transforms into the Einstein Lagrangian 
(3.5), and the KBL gravitational energy-momentum transforms into the Einstein pseu- 
dotensor (3.7). Next, if — > r/^^ and = d'l^^y then the KBL current J('^){0 

goes to i^f^ + Tif in (3.8), the KBL superpotential Jj^^iO transforms into the Freud su- 
perpotential (3.9). That is generally Eq. (3.55) transforms into Eq. (3.10). However, both 
the current and superpotential in Eq. (3.55) are not a simple direct covariantization of the 
quantities (3.7) and (3.9). Indeed, first, the KBL expressions hold on arbitrary curved back- 
grounds, not only on flat backgrounds in curved coordinates. Second, J^'^^ includes the spin 
term (3.51) analogous to the quantity (3.26). It plays its crucial role because then rotation 
Killing vectors can be used giving a reasonable definition of angular momentum. Thus the 
KBL current gives the right ratio of mass to angular momentum for rotating black hole 
solutions. 

The KBL approach received a significant development in applications. Thus, in [178, 
179] the problem of localization was considered. In [180], the conserved quantities and their 
fluxes at null infinity for an isolated sysytem were studied in detail. In the works [175, 181] 
perturbations on FRW backgrounds and conservation laws for them were examined. At 
last, recently [182] the approach was developed for constructing the conserved charges in 
Gauss-Bonnet 5-dimensional cosmology. 

The KBL quantities were also checked from the point of view of the problem of unique- 
ness. Thus, Julia and Silva [183, 184], and independently Chen and Nester [94], stated that 
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the KBL quantities are uniquely defined as associated with the Dirichlet boundary condi- 
tions, under which the action with the Lagrangian (3.41) is variated. This gives evident 
advantages. Both the Lagrangian (3.41) and the energy-momentum (3.50) have only the 
first derivatives. Thus, the Cauchy problem is stated simply. 

Keeping in mind successes presented by the KBL approach it is necessary to note the 
next generic problem related to the canonical derivation. If one chooses different boundary 
conditions, adding different divergences to the Lagrangian, different expressions both for 
currents and for superpotentials appear by the Noether procedure. On the one hand, such a 
situation could be useful and interesting in gravitational theory, when analogies, e.g., with 
thermodynamics are carried out [95]. On the other hand, it is evident that unique expres- 
sions independent on boundary consditions are necessary. As examples, the symmetric 
energy-momentum tensor in classical electrodynamics is such a quantity. The expressions 
in the field-theoretical formulation of GR are also independent on boundary conditions (see 
section 2). In the next subsection we just construct conserved currents and superpotentials 
in the framework of the field approach. 

3.4 Conserved quantities in the field-theoretical formulation 

Already the form of Eq. (3.16) indicates that there is a possible to construct currents and 
superpotentials in the framework of the field approach. Recall that for backgrounds pre- 
sented by an Einstein space the conservation law (2.67) has a place. Then if the Ein- 
stein space has Killing vectors A" one can construct the conserved currents by the same 
way as in the usual field theory (3.23). At the same time the left hand sides both of Eq. 

(2.66) and of Eq. (2.71) contracted with A'^ are presented as divergences of superpotentials 
[68, 100, 144]. However, in the general case of a curved background the conservation law 

(2.67) has no a place. The reason is that there is no an identical conservation both in Eq. 
(2.24): Dy (O^f + d^'^) / and in Eq. (2.32): D^G^'' 7^ 0. The physical reason is 
that the system (2.19) interacts with a complicated background geometry defined by the 
background matter fields $ . 

To construct conserved quantities in the general case we [118] combine the technique 
of previous subsection with the results [100]. Therefore, in the framework of the field 
approach, we are going to the KBL independent variables of the bimetrical approach by 
changing l^^^ g^'^ — g^^ . Then the dynamical Lagrangian (2.22) transforms into 
2kCg2- After that we connect it with the KBL Lagrangian Cq (3.41): 

Cg2 ^Cg- Cgi ^ ^r'' {^%^% - (3.56) 

where Cgi = —{2K)~^{g^'^ — g^'^)R^i,. Returning to the comparison with the Rosen 
Lagrangian Cr [136] we note that jC,g2 is a direct generalization of £r to arbitrary back- 
grounds, whereas jC-g is reduced to jCr for the Ricci-flat backgrounds. 
Now, consider the identity 



£^jC.G2 + d^ieCG2)^0 



(3.57) 
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for the Lagrangian (3.56) and transform it identically to 

dCG2 ^ 



pa- 



pa 



pa 



The quantity 



2^ _ 



dCG2 ,u , d^G2 _ 

gp<y\x + —n — - g, 



pa 



(3.58) 



(3.59) 



is antisymmetric in /z and ly. Thus on the right hand side of (3.58) the quantity [—"^S^^^^] 
plays the role of a superpotential. Besides, the expression (3.59) is connected with the spin 
term (3.51) as 

One has also to note that in (3.58) 



5Cg2 _ 

5g 



pa 



(3.61) 



it is exactly the expression defined in (2.25) and (2.26). 

Now, substitute jCg2 = — -^Gi into Eq. (3.58) and subtract it from the identity (3.53) 
where j'* is defined in (3.45). Returning to the variables l^'^: g^'' — g^^ — > l^^ , one obtains 
the new identity: 

^G^^r + -f^Rxvi" + = Dj^^ (3.62) 
with the new superpotential connected with the KBL superpotential J^s^ given in (3.54) as 



IJ.V _ jfif , 2c IJ-^ cA 
{s) = -^(c) + -^A ? 



(3.63) 

= -iPl^'DpC''^ +V^'%e (3.64) 

K 

= ^ (jP^i'Dp^''^ + (^^'dJ"^" - D^n^'X") ■ (3.65) 

For 5^^, r]f^i, and ^ = S'^^^ the superpotential J^^^ transforms into the the Papa- 
petrou superpotential (3.15) (see the form (3.64)). The last term on the left hand side of Eq. 
(3.62) is connected with the expression (3.52) as 



dCc 



(3.66) 



The way of constructing (3.56) - (3.62) is natural, but a cumbersome one. Formally, 
one concludes that the final identity (3.62) is a result of a difference between identities 
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(3.43) and (3.57). Thus, the Noether procedure could be applied directly to the difference 
Cq — = J^Gi- Indeed, recalculating the identity 

£^CGi+d^{etGi) = ^ (3.67) 

one obtains directly the identity (3.62). The result is not degenerated because Cqi contains 
derivatives of the background metric up to the second order. 

After substituting the field equations (2.24) into the identity (3.62) one obtains the con- 
servation law in the form: 



where the conserved current J^'^^ is expressed through the superpotential J^^^y The general- 
ized total energy-momentum tensor density is 

%). = 4*"*^'' + ^ (^"^^^A. - ^^^) ■ (3.69) 

The role of interaction with the background is played by the term {k)~^ {l^^ R\i, — ^^^). 
Next, if one uses the form of the field equations (2.32), then (3.69) looks as 

%)u = i^u^^^^ + \i^''Rxi^, (3-70) 

and one needs to consider i^^^^^ as the energy-momentum of perturbations and 
{k)~^I^^R\i, as the interacting term. It is interesting to note that the form of the energy- 
momentum (3.70) is more close to the form of the KBL energy-momentum (3.49). The 
other point is that the general term ~ l^^Rx^ destructs the symmetry of ^^gy 

Let us present also gauge invariance properties of the presented current. Substituting 
(2.49) and (2.50) into (3.68) and using the identity (3.62) one obtains for = A": 



1^ dg^" k\ 



k=l 



(3.71) 



thus J^^^ is gauge invariant up to a divergence on the dynamic and background equations. 

Conserved quantities in Eq. (3.68), like the KBL quantities, satisfy all the requirements 
of subsection 3.2. Besides, the field-theoretical derivation has some advantages. First, 
the quantities in Eq. (3.68) do not depend on divergences in the dynamical Lagrangian. 
Second, the current in Eq. (3.68) is defined by the unique quantity, the generaUzed energy- 
momentum without a spin term. Third, the conservation law (3.68) explicitly de- 
scribes perturbations. 

However, the currents and superpotentials in the framework of the field approach bring 
in themselves the Boulware-Deser ambiguity defined for the total energy-momentum tensor 
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in subsection 2.8. Let us outline this problem in detail. Basing on the Lagrangian (2.78) 
one can write out the identity: 



Substituting the field equations (2.80) one obtains 



(3.72) 



(sa) 



where is defined in (3.66) with exchanging l^^'^ by l'^^-^ (see the definition (2.81)). The 



(a) 



generahzed superpotential is 

(sa) ^ V («) ^ W W'^^ J 



(3.74) 



Thus, taking into account the difference in perturbations (2.82) one obtains an analog of the 
Boulware-Deser ambiguity in the definition of the superpotentials: 



In the next subsection we present arguments to resolve this ambiguity. 



3.5 The Belinfante procedure on curved backgrounds 

Returning to the problems of the canonical approach (see the end of subsection 3.3) we re- 
mark that the classical Belinfante method resolves analogous problems in a canonical field 
theory in Minkowski space. Indeed, the Belinfante corrected current (3.29), unlike the cur- 
rent (3.25), is presented without a spin term. At the same time the BeUnfante symmetrized 
energy-momentum is equal to the symmetrical (metric) energy momentum (3.30), and thus 
it does not depend on divergences in Lagrangian. Of course, a perturbed system on a curved 
background in GR is a more compUcated case. However, the KBL model [176] looks an 
appropriate one for an apphcation of the Belinfante procedure. Indeed, for a Killing vector 
the current J^'^^ in (3.48) has the form (3.25) with the correspondent energy-momentum 

0|^^^^ and spin term a^^P'^. Thus, it is anticipated that the Belinfante procedure could 1) 
resolve the problems of the KBL approach, 2) connect it with the field-theoretical approach 
and, consequently, 3) resolve problems of the last also. Therefore, in this subsection, we 
generalize the Belinfante method for constructing conservation laws in the perturbed GR 
on arbitrary curved backgrounds and following the papers [115, 116] apply it to the KBL 
model. 

Thus, for the spin tensor density (3.51) with the rules (3.27) we construct the quantity: 

gfMiyp ^ _gytxp ^ ^p[ai!^] _|_ ^iJ-ipi'] _ ^^[ptA (3.76) 

and present the KBL conservation law (3.53) in the equivalent form: 

J{c) + {S^^'ip) = du {J^^) + S^^'ip) . 0.11) 
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Defining the left hand side as a symmetrized current J^'^^ and the right hand side as a 
divergence of a superpotential J^^^y we rewrite (3.77) in the form: 

•^(B) = ©(6)1/^" + ^{B) = ^^J{By (3.78) 

Here, the spin term is absent; if is a ICilling vector of the background Z-term disappears; 
the symmetrized energy-momentum tensor density has the form: 

^U^ = Ky + DpS>^''.. (3.79) 

Now, consider properties of Belinfante symmetrized quantities. The symmetrized 
energy-momentum tensor density (3.79) has the structure 

Gf^) = i^^ + {f^^t^P - T"^) + l-iP-R^,g^^- + . (3.80) 

Here, the first term is the symmetrical energy-momentum tensor density for the free gravi- 
tational field: 

+ g^^ (a-,a('; - A^^A(^, - A^^A(^,) (3.81) 

The second term in (3.80) is the perturbation of the matter energy-momentum, the third and 
fourth terms describe the interaction with the background geometry. 

Because the fourth term in (3.80) is antisymmetric, the energy-momentum is symmetric: 
(B) ~ {B) ^'^^y ^f^'^ ~ ^d^Li^' ^'^^ cases when backgrounds are Einstein 

spaces in Petrov's terminology [143]. Thus, the Belinfante symmetrization is restricted in 
the sense of constructing really symmetric quantities. Similary, the differential conservation 
law Di,Q'^^^ = has also a place if and only if backgrounds are Einstein spaces. However, 
even on arbitrary curved backgrounds for Killing vectors the current is defined, like 
in (3.29): J^^^^ = Q^^By^"" , and is conserved: {p^(B)u^'') = 0- Thus both "defects" 
compensate one another. This could be useful, e.g., for constructing angular momenta of 
relativistic astrophysical objects on a FRW background, which is not an Einstein space, but 
which has rotating Killing vectors. 

Unlike the KBL quantities, second derivatives of g^^y appear in the energy-momentum 
tensor density (3.81). This needs some comments also. The canonical energy-momentum 
i^, see (3.50), is quadratic in first order derivatives, and this is the normal behaviour of 
a conserved quantity relating to initial conditions. Consider the local quantities J^'^^ in 
(3.77) and (3.78). Suppose that initial conditions are defined on a hypersurface at a given 
time coordinate = const. Then it is important to examine initial conditions for zero's 
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component J(°j) = Jj^ + dk{S^'"'^a) only, see Eq. (3.33). Such a form is because Si^"" 

is anti-symmetric in the first two indices, see (3.76). Thus since J^^^ and S'^'^^^a contain 

only first order time derivatives, J^^^ itself contains only first order time derivatives of the 
metric and therefore has the normal behaviour with respect to initial conditions. Thus, this 
requirement may be unnecessarily restrictive. 

Now let us discuss the properties of the quantities in (3.78). Comparing the expressions 
in Eqs. (3.60) and (3.76) and definitions in Eqs. (3.63) and (3.77) we find out that the 
superpotentials are identical. The same conclusion is related to Z-terms, thus 

= ig, (3.82) 

= ^r^)- (3-83) 

Next, taking into account these relations and comparing (3.68) with (3.78) we conclude that 
"^(B) ~ "^{s) consequently it has to be 

©fs). = ■ (3-84) 

Although the structure of these generaUzed energy-momenta are different, direct calcula- 
tions with the use of the field equations also assert the claim (3.84). Thus, on arbitrary 
curved background one cannot use separately the energy-momentum of the free gravita- 
tional field neither t^j,, nor t^^. One needs to use the total energy-momentum tensor den- 
sity, which is the same, as is seen from in (3.84). 
From all of these one concludes 

• The construction of the conservation laws for perturbations on arbitrary curved back- 
grounds in GR with the use of the generalized Belinfante symmetrization transforms 
the conservation laws of the canonical system into the ones of the field-theoretical 
approach. Thus, the Belinfante method is a "bridge" between the canonical Noether 
and the field-theoretical approaches. 

In spite of that Eq. (3.84) is an analog of the relation (3.30) in a simple field theory (3.18), 
this conclusion is not so evident. Indeed, we cannot apply the Belinfante method only 
in the framework of the KBL bimetric system (3.40) because we cannot define a sym- 
metrical energy-momentum analogously to (3. 19). Really, the variation of the Lagrangian 
(3.40) with respect to the background metric gives only the background Einstein equations, 
whereas the variation of the Lagrangian (3.41) gives only the non-dynamical background 
Einstein tensor. Only a non-trivial connection with the field approach turns out fruitful 
leading to (3.82) - (3.84). 

Relations (3.82) - (3.84) solve the aforementioned problems of the canonical construc- 
tion. Indeed the quantities 0(s)i/' ^{b) ^{b) ™ quantities in the field- 
theoretical frame and therefore they do not depend on divergences in the Lagrangian. A 
direct mathematical substantiation is given in subsection 4. 1 . Thus we have arrived a one 
of desirable results. 

Now, return to the ambiguity of the field approach in currents and superpotentials noted 

at the end of subsection 3.4. Remark that the KBL approach and the Belinfante symmetriza- 
tion do not depend on a choice of the variables, hke g^^ , g^^, ^/--gg^^ , ■ ■ ■, and uniquely 
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lead to the conservation law (3.78). On the other hand, it coincides uniquely with conser- 
vation law (3.68) defined by the second decomposition in (2.73). Thus it is a theoretical 
argument in favor of a choice I'^^s^ = l'^'^ in (2.81). To support this theoretical conclusion 
below we present a useful test. Recall that differences between superpotentials in the family 
(3.84) appear in the second order. Just the second order is crucial, e.g., in calculations of 
energy and its flux [185] at null infinity. As it was checked in [116], only the superpotential 
(3.65) from the family (3.75) gives the standard result [185], whereas all the others, e.g., 
the superpotential corresponding to the first decomposition in (2.73) and defined in [68], do 
not. 

3.6 Applications 

The above presented formaUsm, its technique and expressions (in particular, the energy- 
momentum tensor and its properties, including gauge invariance ones) were used in 
various appUcations in GR. Shortly we note that tlf"*"* was used in development of quantum 
mechanics with non-classical gravitational self-interaction [108, 109], in the framework of 
which inflation scenario was analyzed [110, 111]. Some other applications are reviewed 
below in more detail. 

3.6.1 The weakest falloff at spatial infinity 

The proof of the positive-energy theorem in [70] - [73] has renewed a great interest to 
asymptotically flat spacetimes in GR. In particular, these investigations showed that the 
standard asymptotic conditions at spatial infinity (1/r falloff in metric) can be significantly 
relaxed (see, e.g., [186] - [189] and references there in). These models are very appropriate 
for applications in the framework of the field approach. Indeed, asymptotically a back- 
ground is just chosen naturally, it is a Minkowski space at spatial infinity. Gravitational 
field at spatial infinity is weak metrical perturbations with respect to this background. We 
have done such applications: both inLagrangian [112] and in Hamiltonian [113] derivation 
we examine the weakest asymptotic behaviour at spatial infinity. We consider so-called real 
isolated systems, for which all the physical fields are effectively concentrated in a confined 
space at finite time intervals, and a falloff of which is invariant under asymptotic Poincare 
transformations. 

As an initial falloff we use the standard one: 1/r. Then choosing a background space- 
time as the Minkowski space we evaluate a field configuration in the framework of the field 
approach. For simplicity it is assumed that a manifold, which globally supports the physical 
metric, supports also the auxiliary flat metric. To search for the weakest falloff global con- 
served quantities of the isolated system (integrals of motion) are examined. We construct 
them in the form of the surface integrals (3.39) with the use of the superpotential defined in 
(3.65) and the ten Killing vectors of the Minkowski space. 

To find the weakest falloff conditions one uses the gauge invariance properties of the 
field formulation of GR in section 2.6. Taking into account the general formula (3.71) for 
the current (3.68) we conclude that the integrals of motion are gauge invariant up to surface 
term on the solutions of the field equations. Then we search for the weakest falloff for the 
gauge potentials and their derivatives which ensure vanishing the non-invariant terms. 
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After that we find out the weakest falloff for the gravitational potentials: 

11^" = 0+{r-^) + 0-{r-^); e + S > 2 , 1 > e > ^ , S > 1 (3.85) 

where (+) and (— ) mean even and odd functions with respect to changing sign of the 3- 
vector n = x^/r. To obtain the simple conclusion that the falloff can be weakened up to 
r"^/^ one needs to examine energy and momentum only (see, e.g., [187, 188]). To obtain 
the more detailed falloff (3.85) one has to study all the ten integrals of motion. In this re- 
lation we remark also the works [186, 189]. In [186] the algebra of asymptotic Poincare 
generators is studied, and in [189] the finiteness of the ADM definition of angular momen- 
tum is examined. The result (3.85) based on the gauge invariance requirement coincides 
with the results in [186] corrected not significantly, and with the results in [189], if the 
last are adopted by some natural requirements (see [112, 113]). Thus the three different 
approaches give the same weakest falloff. 

Historically, among the ten integrals of motion, not enough attention was payed to the 
Lorentzian momentum. In this relation the paper by Regge and Teitelboim [190] who in- 
cluded the falloff of odd and even parity, and the paper by Beig and OMurchadha [191] who 
improved the results of the work [190] were the most important earlier works. Currently 
this gap is closing [117], [192] - [194]. In [117] we demonstrate that the quasi-local Brown- 
York [79] center of mass integral asymptotically agrees with the one given in [191] and with 
the one expressed by the curvature integrals, e.g., in [195]. In [192] the falloff of the grav- 
itational potentials is analyzed in n -I- 1 dimensions with special attention to the Lorenzian 
momentum and matter variables. In [117] and [192] (for n = 3) the falloff exploited by 
Beig and OMurchadha [191] is used. It is r^^ for the even parity first term and r^^, 6 > 1, 
for the next of arbitrary parity term. Comparing with the spectrum of the conditions (3.85) 
one finds that this condition corresponds to £ = 1 and (5 > 1 in the full spectrum (3.85). As 
an example, the other end of the spectrum is e > 1/2 and 6 > 3/2. In [193] and [194], in 
the framework of the covariant Hamiltonian approach [94] the quasi-local center of mass 
integral in the teleparallel gravity and GR is considered. This approach gives quasi-local 4- 
momentum, angular momentum and center of mass integral in the framework of an unique 
relativistic invariant (Poincare invariant) description. 

3.6.2 Closed worlds as gravitational fields 

In the paper [101], we show how a closed- world geometry can be viewed as a gravitational 
field configuration, superposed on a topologically trivial geometrically flat spacetime. As 
seems, it is not very natural. However, we do this without contradicdons. The space compo- 
nents of this configuration corresponds to a so-called stereographic projection of 3-sphere 
onto a flat 3-space E^, when the "south" pole corresponds to the origin of the coordinate 
system, whereas the "north" pole is identified with all the points at infinity of E^. Thus, to 
make the full identification the "north" pole is "knocked out" of the sphere to simplify 
its topology. 

The background Minkowski space has an auxiliary character, as it has to be. Formally, 

the field configuration is specified in an infinite volume. However, with the use of the phys- 
ically reasonable measurements examining hght signals in gravitational fields one finds the 
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standard volume for (see [1]). By considering the metric relations established through 
the physical measurements an observer will infer that homogeneity, an equivalence of ob- 
servers at all points of space, has a place. 

A picture when topological and geometrical properties are replaced by the properties of 
an effective field propagated in a trivial spacetime could be interesting and useful. Thus, 
Rubakov and Shaposhnikov [40] have shown that not too energetic scalar particles could 
become effectively trapped in a potential well even in a topologically trivial universe, al- 
though nontrivial classical solutions would have to be present to play a role of an external 



The field configuration presenting the closed world is static and has an asymptotic be- 
haviour This ensures zero energy, momentum and angular momentum of the system, 
which then could be treated as a microuniverse. Indeed, such characteristics are true for 
a Minkowski vacuum where classical fields and particles are absent, and thus an idea of 
quantum birth of the universe could be supported. 

3.6.3 A point particle in GR 

The Schwarzschild solution is the one of the most popular models in GR. Frequently, in 
spite of a complicated geometry, it is treated as a point mass solution in GR [1]; currently 
the interest to this viewpoint is renewed [196] - [199]. However, if one considers GR in the 
usual geometrical description, then this interpretation meets conceptual difficulties (see, for 
example, [200] and also [114]). 

Except a pure theoretical interest the point particle description could be interesting and 
useful for experimental gravity problems. Gravitational wave detectors such as LIGO and 
VIRGO will definitely discover gravitational waves from coalescing binary systems com- 
prising of compact relativistic objects. Therefore it is necessary to derive equations of mo- 
tion of such components, e.g., two black holes. As a rule, at an initial step the black holes 
are modeled by point-like particles presented by Dirac's 5-function. Then consequent post- 
Newtonian approximations are used (see, for example, [201, 202] and references therein). 

The aforementioned difficulties [200], at least, do not appear in Newtonian gravity. 
To describe a point particle one has to assume that the mass distribution has the form 
p{r) = m6{r) where (5-function satisfies the ordinary Poisson equation, which in spher- 
ical coordinates is 



Then, the Newtonian potential cj) = m/r satisfies the Newtonian equation in the whole 
space including the point r = 0. Besides, the whole mass of the system is defined by 



particle density can be integrated providing the standard mass m. 

In [114]) we show that, analogously to the Newtonian prescription, the point mass in 
GR can be described in a non-contradictory manner in the framework of the field formu- 
lation. The Schwarzschild solution was presented as a gravitational field configuration in 
a background Minkowski space described in the standard spherical Schwarzschild (static) 
coordinates. The concept of Minkowski space was extended from spatial infinity (frame of 



field. 




(3.86) 



the unique expression fj, dx^p{r) 



where both the usual regular distribution and the point 
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reference of a distant observer) up to the horizon r = rg, and even under the horizon in- 
cluding the worldline r = of the true singularity. The configuration satisfies the Einstein 
equations at all the points of the Minkowski space, including r = 0, if the operations with 
the generalized functions [203] are valid. The energy-momentum tensor was constructed. 
Its components are presented by expressions proportional to 5(r) and by free gravitational 
field outside r = 0. The picture is clearly interpreted as a point particle distribution in GR. 
Indeed, the configuration is essentially presented by (5-function, one can use the volume 
integration over the whole Minkowski space and obtain the total energy mc? in the natural 
way. However at r = rg the energy density and other characteristics have discontinuities. 
A"visible" boundary between the regions outside and inside the horizon exists and does not 
allow to consider an evolution of events continuously. 

However, it is not a real singularity. In the field formulation this situation is interpreted 
as a "bad" fixing of gauge freedom, which can be improved. That is the break at r = can 
be countered with the use of an appropriate choice of a flat background, which is determined 
by related coordinates for the Schwarzschild solution. At least, the use of the coordinates 
without singularities at the horizon could resolve the problem locally at neighborhood of 
r = Tg. Next, it would be natural to describe the true singularity by the world fine r = 
of the chosen polar coordinates. Besides, it is desirable to have appropriate coordinates, in 
which the Schwarzschild solution conserves the form of an asymptotically flat spacetime. 

Recently Pitts and Schieve [130] defined and studied properties of a so-called "77- 
causality". Its fulfilment means that the physical light cone is inside the flat Ught cone 
at all the points of the Minkowski space. It is necessary to avoid interpretation difficulties 
under the field-theoretical presentation of GR. By this requirement all the causally con- 
nected events in the physical spacetime are described by the right causal structure of the 
Minkowski space. A related position of the light cones is not gauge invariant. We consider 
this requirement only to construct a more convenient in applications and interpretation field 
configuration for the Schwarzschild solution. The requirement of the ry-causality can be 
strengthened by the requirement of a "stable ry-causality" [130]. The last means that the 
physical light cone has to be strictly inside the flat Ught cone, and this is important when 
quantization problems are under consideration. Indeed, in the case of tangency a field is 
on the verge of ?7-causality violation. Returning to the presentation in the Schwarzschild 
coordinates in [114] we note that it does not satisfy the ?7-causality requirement. 

In the paper [120], taking into account the above requirements we have found a desir- 
able description. A more appropriate gauge fixing corresponds to the stationary (not static) 
coordinates presented in [105, 106] and recently improved in [130]. We consider also the 
contracting Eddingtom-Finkenstein coordinates in stationary form [3]. These two coordi- 
nate systems belong to a parameterized family where all of systems satisfies all the above 
requirements. The transformation 



gives just the parameterized by a G [0, 2] family of metrics. The cases a = 1 and a = 
correspond to the first and to the second examples above. The stable ??-causaUty is not 
satisfied with a = 0at0<r<oo. Thus, all the configurations a € (0, 2] are appropriate 
for the study both classical and quantum problems, whereas the case a = could not be 




(3.87) 
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Figure 4: The energy distribution for the field configuration associated with the 
Schwarzschild solution. The case: a G (0, 2]. 

useful for the study quantized fields. Properties of field configurations corresponding to 
a G (0, 2] qualitatively are the same as for a = 1. In the terms of the field approach [100], 
all the field configurations for a G [0, 2] are connected by gauge transformations and are 
physically equivalent. 

Both the gravitational potentials and the components of the total energy-momentum 
tensor for the field configurations with a G [0,2] have no a break at r = Vg. The energy 
distribution is described by the 00-component of the energy-momentum tensor, and for 
a G (0, 2] qualitatively is presented on the figure 4. Then the total energy of the system 
can be calculated both by the volume integration and by the surface integration over the 
2-sphere with r ^ oo giving m(?. For the case a = the components of the total energy- 
momentum tensor for the field configuration have the simplest form: 

4°o* = mc25(r), 
4f = -mc25(r), 

^'ab = -^VABmc^Sir). (3.88) 

This energy-momentum is concentrated only at r = 0, see the energy distribution on the 
figure 5. The other component t\°i and the angular ones formally could be interpreted 
as related to the "inner" properties of the point. Indeed, they are proportional only to 6{r) 
and, thus, describe the point "inner radial" and "inner tangent" pressure. At last, again the 
total energy for (3.88) is mc^. 

The field configurations for a G [0, 2] satisfy the Einstein equations (2.66) at all the 
points of the Minkowski space including r = 0. Then, keeping in mind the presentations 
on the figures 4 and 5 one can conclude that an appropriate description of a point particle in 
GR is approached. It is directly and simply continues the Newtonian derivation and, at the 
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Figure 5: The energy distribution for the field configuration associated with the 
Schwarzschild solution. The case: a = 0. 



same time, does not contradict GR corresponding its principles. 
3.6.4 Integral constraints in FRW models 

In this subsection we consider the linear perturbations on FRW backgrounds from the point 
view of the conservation laws (3.78), or the same (3.68), for the detail calculations see [116]. 
Write the background metric ds"^ in dimensionless coordinates = r],x'' as follows [153]: 

ds"^ = g^^dxf'dx'' = a^{'q){drf - ftidx'^dx^) = a'^{r])e^udx^'dx'' , (3.89) 

a(r?) is the scale factor, fki = 6ki + Mfcm(^i„x'"x"'(l - kr'^)~^, f = det fki = (1 - fer^j^^ 
/c = or ±1, and = 6kix^xK We define the perturbations by gfj,i, = c^'iey.v + ^ij.u) 
as often used in cosmology (see, e.g., [126, 204]), which are connected with the the usual 
definition (2.13) in linear approximation: 

= c?^f[-e^Pe''' + ie^'^e^") V • (3-90) 

Because the FRW spacetimes are conformal to Minkowski space the conformal Killing 
vectors of FRW spacetimes are thus those of Minkowski space written in t]^ x^ coordinates 
(see [205] - [207]). There are 15 conformal Killing vectors, ^a, ^ = 1, 2, . . . , 15, some 
of which are pure Killing ones. 

These vectors we use as displacement vectors in the expression (3.78). Next, we make 
linear combinations of ^yi's with functions of rj, say c^{ri)^A- In general the C^^a's are not 
conformal Killing vectors but some linear combinations turn out to be rather simple with 
clear physical interpretations in appropriate gauge conditions. Integrating the conservation 
law (3.78) we construct the integral relations connecting perturbations inside of a volume 
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with perturbations at its surface. Thus we consider a sphere r = const at constant time rj 
and integrate zero components of (3.78) J^^^ = diJ^^^ over such a ball. After cumbersome 
calculations for perturbations (3.90) one has 15 volume integrals connected with surface 
integrals. 

Besides the matter perturbations ST^ the volume integrands contain two field quantities 
Q = 5{—Dnn^) — the perturbation of the external curvature scalar of 77 = const, and 
— the space trace of the metric petrurbations. There are four linear combinations which do 
not involve /i™; they are associated with some linear combinations of the conformal Killing 
vectors, which are just Traschen's [172] "integral constraint vectors". 

Of course, we must still impose four gauge conditions to fix the mapping of the per- 
turbed spacetime onto the background. One gauge condition that simplifies almost all 
volume integrands is the "uniform Hubble expansion" gauge Q = discussed, e.g., by 
Bardeen [10]. Then, 14 of the 15 volume integrands reduce to Unear combinations of 5Tjj 
only, and thus the 14 relations present "integral constraints", among which the 4 are new. 
The integrals represent momenta of the matter energy-momentum tensor of order 0, 1 and 
2 in powers of when A; = and with similar interpretations when k = ±1. There re- 
mains one integral that contains both ST^ and h!^: it is related to so-called conformal time 
translations for A; = ±1 or to so-called time accelerations for A; = 0. 

Another gauge condition often used (see, e.g., [204]) is Vj^j, = in which hi. is the 

traceless part of /i'^ ; these gauges remove the longitudinal modes of the gravitational waves. 
Combining V/ZiJc = with Q = one finds four relations that are independent of the 
gravitational radiation. 

4 D-dimensional metric theories of gravity 

As was discussed in Introduction, multidimensional gravitational models become more and 
more popular, their solutions (for example, generalized p-branes and brane-world black 
hole solutions [208] - [211]) induce an arising interest. Other generalizations, say, scalar- 
tensor theories of gravity [212], are developed intensively also. In this section, keeping 
in mind such generalizations we construct a perturbed scheme and conservation laws for 
perturbations in a generic metric theory. The presentation is based on the works [103, 104, 
119, 121]. Besides, here the approach is presented in the united scheme, and more details 
are given. 

4.1 The main identities 

In this subsection, we present necessary identities and stress their important properties. Let 
the system of fields, set of tensor densities Q^, be described by the Lagrangian 

£ = £(Q^;Q^,a;Q^,„;3) (4.1) 

including derivatives up to the second order. We assume that it is an arbitrary metric theory 
of gravity in D dimensions with the generahzed gravitational variables {g"',"^^) and the 
matter sources . Thus = {(5", ^^), ^'~^}; the variables are included keeping 
in mind, say, scalar-tensor or vector-tensor theories; the metric variables g"- are thought 
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as defined in (2.74), only in D dimensions. Then we include as usual. The ordinary 

derivatives da are rewritten over the covariant ones by changing Qb,t = D-j-Qb — 
Qsla- Then the Lagrangian (4.1) takes an exphcitly covariant form: 



(4.2) 



After that indexes are shifted by g^^ and g^^". 

Now we use the standard technique (see, for example, [140]). For the Lagrangian (4.2), 
as a scalar density, we write out again the identity: £^jC,c + (^">Cc),a = 0, which can be 
rewritten in the form: 



sqb "^"^ ysQB 

In (4.3), the coefficients are defined by the Lagrangian in unique way: 



0. 



(4.3) 



Da - 



DaQs 



diD^DaQs 



(4.4) 



m„ = 



(- 



Qb\1 



d{DrDaQB) d{DpDaQB) 



diDf^DaQB) d{DrDaQB) 



(4.5) 



(4.6) 



The coefficient at ^'^ in the first term in (4.3) is identically equal to zero (generalized 
Bianchi identity). Thus the identity (4.3) has the form of the differential conservation law: 



DJ'' = doT = 

with a generalized current 

r = - [u^r + K'^Dri" + K^^DpDrC 

We also use the helpful form: 

r ^ - [{u- + hf^R\^,)e + mfr'dyi^ + 
with z-term defined as 



(4.7) 



(4.8) 



(4.9) 
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Again, if = A", then = and the current (4.9) is determined by the energy- 
momentum {u + ni?) -term and the spin m-term. Opening the identity (4.7) and, since 
da^"^, dpcti"^ and d^paC" arbitrary at every world point, equating independently to zero 
the coefficients at Da^", D(^pa)i" and -^^(7/3a)€'^ we get a "cascade" (in terminology by 
JuUa and Silva [183]) of identities: 

J^a "cr 2 A -"-o- po 3 ""A ^l^a pa — 

4"^'^^ = 0. (4.11) 

The above constructions are the generalization to arbitrary curved backgrounds of the ex- 
pressions given by Mitzkevich [140]. 

Since Eq. (4.7) is identically satisfied, the current (4.9) must be a divergence of a 
superpotential z"'^, for which dj^aV^^ = 0, that is 



= Dpi"^ = d^r^. (4.12) 



Indeed, substituting vf^ and rh'^^ from Eqs. (4.11) directly into the current and using alge- 
braic properties of n"^''^ and R'^ppfj, and the third identity in Eqs. (4. 1 1) we reconstruct (4.9) 
into the form (4.12) where the superpotential is 

= (pAni"'^!^ - mi"'^!) T - In^^^^^Dx^" ■ (4.13) 

It is explicitly antisymmetric in a and /3. Of course, Eq. (4.12) has also a sense of the 
differential conservation law Eq. (4.7). 

Let us find contributions into currents and superpotentials from a divergence in the 
Lagrangian 5dCc = div = (F ^y. For the scalar density dy ^i, one has the identity {£^(P + 
i'^dy ,u).a = 0, which gives contributions 5^1°^ into the current (4.8) and 5^?"^ into the 
superpotential (4.13). The additional quantities are 

<5d< = 2:D^(4"d^l), (4.14) 

5dmf = 24"d^l, (4.15) 

5dhf^ = 0, (4.16) 

5di"'^ = -2^^"d^^. (4.17) 

Then, for £-c + ^d^^c Eq. (4.12) changes as 

r + Sdi" = Dp (f^ + jrf^^) , (4. 18) 

where the changes do not depend on a structure of d^. 

Next, using the general Belinfante rule (3.27) we define a tensor density 
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add D p{s°'l^^ ^cr) to both sides of (4.12), and obtain a new identity: 

rB^D^e/^df^if. (4.20) 
This modification cancels the spin term from the current (4.9): 

i% = (-< - nf ^i?'^^, + D^r^) r + m) ^ i^iaC + m), (4.21) 

and a new z-term disappears also on KilUng vectors of the background: 

Z%{0 = (y'^^f + ^Tf' - ^l^r') Cr/3 + (2D^0Cr) - DaQrg''') • (4.22) 

Thus, the current i% is defined, in fact, by the modified energy-momentum tensor density 
u%^. Because the new superpotential depends on the n-coefficients only: 

if ^ 2 (ppnl-'^lP + DrhY^''-g^^%,) ^ - Ih^^^^'D,^^, (4.23) 

then due to the definition (4.6) it vanishes for Lagrangians with only the first order deriva- 
tives. On the other hand, the superpotential (4.23) is well adapted to theories with second 
derivatives in Lagrangians, like GR or the Einstein-Gauss-Bonnet gravity. It is also impor- 
tant to note that for the superpotentials (4.13) and (4.23) their forms do not depend explicitly 
on a dimension D; the same property is related to the Deser and Tekin superpotentials [75]. 

It is important to note that the Belinfante procedure applied to (4.18) cancels the quanti- 
ties Sdi" and 5^?"'^ and gives again (4.20). Let us show this. The quantity (4. 19) constructed 
for Sdrh^^ in (4.15) gives 

Sds'^'^'^^a = 2^[°d^l . (4.24) 

Then, adding Dp{s°'^'^^a + Sds"'^'^ ^a) to Eq. (4.18) one cancels completely the spin term 

and suppresses Sdu"^ and 5dt"^ ■ Indeed, combining (4.14) and (4.17) with (4.24) one has 
5dK - D(^{5dS°'^f'gpa) = and Sdi""^ + Sds'^'^'^^a = for an arbitrary d". This just 
supports the claim that J(^By ^^{b)v and J^'^^ in (3.78) are independent on divergences. 

4.2 The field-theoretical formulation for perturbations 

In this subsection, following the method of section 2 we develop a perturbed derivation. 
The field equations corresponding to the system (4.1) are derived as usual: 

5C/5Q^ = {). (4.25) 

Let us decompose Q onto the background part Q and the dynamic part q : 

Q^ = Q'^ + q'^. (4.26) 
The background fields satisfy the background equations 

51/ 5Q^ = (4.27) 
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where L = C{Q). Taking into account the form of the Lagrangian (2.19) (or (2.78)) we 
will describe the perturbed system by the Lagrangian 

£<^2'"(Q, q) = C{Q + q) - -2 + div (4.28) 

6Q 

instead of the original Lagrangian C{Q) = C{Q + q). As before, the background equations 
should not be taken into account before variation of C'^y^{Q, q) with respect to Q^- Using 
the evident property 6C{Q + q)/5Q = 5C{Q + q)/Sq^, the equations of motion related 
to the Lagrangian (4.28) are presented as 

-r^ = ^ jC{Q + q)-jC =0. (4.29) 

Sq^ SQ^ L J 

It is clear that they are equivalent to the equations (4.25) if the background equations (4.27) 
hold. 

Defining the "background current" 



5Q 5Q 5Q 



and combining this expression with (4.29) one obtains another form for the equation (4.29): 

5 s 5I 
SQ SQ 



GA' + ^A'^--L<i''^ = t\. (4.31) 



The left hand side in (4.31) is a linear perturbation of the expression SC/SQ^ in (4.25), 
G\ is a pure gravitational part. Really gravitational variables g*^ and ^'^ have a different 
nature. Conservation laws are connected with symmetries of a geometry related directly to 
g"-, not with Thus, we will consider as included in setting thus in calculations 
= and = {g^, $^}. Later we will consider the lagrangian 

C = CDg{g'')+tDra{g\^'') (4.32) 

with the pure metric gravitational part Cdq- Then the equations (4.30) are separated into 
gravitational and matter parts as follows 

$g^--4-g^%^=i?;. (4.34) 

As is seen, the form of the perturbed equations (4.31), (4.33) and (4.34) is a quite universal 
form with the "background current" as a source on the right hand side. Contracting (4.33) 
with 2dg°' /g'^'' one obtains 

^ "I" ^ ixv ~ ^liv (4.35) 
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where the linear operator is defined as 



(4.36) 



with independent gravitational variables defined as in (2.81). The right hand side in 
(4.35) is defined by variation of (4.28) 

=2 ^^^'"^^'^^ (4 37) 

and is the generalized symmetric energy-momentum. 

The equations (4.35) generalize the equations of GR (2.24), they generaUze also the 
Deser-Tekin equations [75, 213, 214] constructed for the quadratic theories by direct cal- 
culations. Below we will show that for the vacuum backgrounds (when = 0) 
D^G^^ = 0, and thus the energy momentum tensor f^^^ is differentially conserved, like 
in the Petrov spaces (2.67) in GR. The presented here model has the same properties for ex- 
pansions as was described in subsection 2.5 and the same gauge properties as was described 
in detail in subsection 2.6. 



4.3 Currents and superpotentials in the field formulation 

In the recent series of the current works [75, 76, 213, 214] Deser with coauthors develop 
a construction of conserved charges for perturbations about vacua in metric quadratic (in 
curvature) gravity theories in D dimensions. They apply the Abbott and Deser procedure 
[68] and develop it. The aim of the present subsection is to suggest an approach, which gen- 
eralizes the Deser and Tekin constructions. We construct conserved currents and superpo- 
tentials corresponding to the equations (4.35) derived for an arbitrary metrical gravitational 
theory, and not only on a vacuum background. Below taking into account the definition 
of the gravitational part of hnear operator (4.36) we demonstrate that conserved quantities 
of the system and their properties can be obtained and described analyzing only the scalar 

density Ci = i'^^^{SjCDg/^g°' ), which is the gravitational part of second term in the La- 

grangian (4.28). As an important case we consider explicitly only such theories where Ci 
has derivatives not higher than of second order, like the Einstein-Gauss-Bonnet gravity. In 
principle, our results can be repeated when Ci has derivatives of higher orders, like in [75]. 

Keeping in mind that £i is the scalar density we again follow the standard technique 
and use the results of the subsection 4.1, which are universal. One transforms the identity 
£^£.1 + daiSf^ti) = into the identity: 



D.t = dj't = (4.38) 



where 
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It is just the current (4.8) for the Lagrangian Ci. 

In the case of a vacuum background one has SCng/Sg^'^ = and Ci = 0. Then 
assuming arbitrary Killing vectors = A" one transforms the identity (4.38) into 

Recalling that for the vacuum case = ^rid taking into account the identity (4.41) 
with the definition (4.36) in the equations (4.35) one gets D^t^'^ = 0. This generalizes the 
results in [213, 75] for the quadratic theories. Recall also that in the Lagrangian (2.69) A- 
term was interpreted as "degenerated" matter, then with this assumption the identity (2.72) 
has been approached. However, the Lagrangian (2.69) fully is, of course, a kind of metric 
Lagrangians, therefore the identity (2.72) also is interpreted in the terms of the general 
formula (4.41). 

Again, because (4.38) is the identity the current has to be presented through a di- 
vergence of a superpotential. To construct the last we use the results of the subsection 4.1 
adopted for the Lagrangian £i. We set g^^^ g^,y and construct the coefficients (4.5) - 

(4.6) with Ci = Ci{Q- d^Q; d^M) = D^^Q; D^M), where = {i^, 

Then we go back, 3^1/ g^^, and obtain simple expressions 

m,^^-^ = 2D, (^^] , nr = '^^^^9,^ ■ (4-42) 
\ogpu,^xJ ^pMx 

Then, substituting these into the expression (4.13) we can define the superpotential i'^gy 

which is evidently linear in I'^K^ . In the case of the Einstein gravity it is the superpotential 
(3.74). Thus, we have the identity 

h — c'i^-i-(s) 

which generahzes the identity (3.72) and the superpotential has the form 

^ (i^An'f ' - mlf ) e - fnifl ^^A^ . (4.44) 

Using in (4.39) the field equations in the form (4.35) and the definition (4.36) we write 
out the physically real current: 1^^^ : 

TZ = %\ r + (4.45) 
where the generaUzed energy-momentum is 

^^J. = fff^ _ \ _ i X^tfP^ ^^Dg I pa ^ 5C,Dg 



it = > (4-43) 



' ^ 6g 

Thus finally the identity (4.43) transforms to the conservation law: 

-is) - ^""-^is) ■ 



(4.46) 



^L=dulZ- (4.47) 
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The expressions (4.45) and (4.46) are based on the symmetrical energy-momentum 
(4.37) and generalize the correspondent expressions in GR: (3.68) and (3.69), therefore 
we choose the subscript "{s)". Thus, this subsection, generalizing the results of the subsec- 
tion 3.4 develop the idea to use the identity (3.67) in GR. The expression (4.47), on the one 
hand, generalizes the conservation law (3.73) in GR, on the other hand, it generalizes the 
Deser-Tekin expressions [75, 213, 214]. In the case of a vacuum background and using the 
KilUng vectors = A" one has for the current 

^{s) = i^i^^" ' (4.48) 
that also is in correspondence with the general definition (3.23) in a field theory. 



4.4 Canonical Noether and Belinfante symmetrized currents and 
superpotentials 

The expressions presented in subsection 4.1 are maximally adopted to construct both 
Noether canonical conserved quantities in the framework of the bimetric formulation and 
Belinfante corrected quantities. To construct such quantities one has to consider a pure met- 
ric part jO-Dg of the Lagrangian (4.32). At the beginning we construct the Noether canonical 
quantities. We construct the corresponding to Cog coefficients (4.4) - (4.6), u'^^, rh'^^, 
hg^'^, with the use of which we present the related identity (4.12): 

= dfsif . (4.49) 

Next, following to the KBL ideology (3.41) we construct a metric Lagrangian for the per- 
turbed system: 

CDG = ^Dg-lDg + dJ''. (4.50) 

Then, we apply the barred procedure to the identity (4.49) and take into account the diver- 
gence keeping in mind (4.14) - (4.18). As a result we obtain the identity corresponding to 
(4.50): 



/9 



(4.51) 



After substituting the dynamical equations (4.25) and the background equations (4.27) 
in the form 

(4.52) 
(4.53) 



Sj^Dg 


S^Dm 






SCog 


SC-Dm 


sr 


sr 



into u^^ and Ug^ one obtains Ug^ and U g^j, respectively. Then the identity (4.51) transforms 
into a real conservation law: 

Jf,)(0=a^X(t(e). (4.54) 
The left hand side (current), in correspondence with (4.9), is 

^ ^."(c)r + S'(!)'d[0^ + (4.55) 
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where the generalized canonical energy-momentum, spin and Z-term are 

^ -[{6U^„ + 2Dp{6^^d^^))+6hf^''R'p^,] (4.56) 

S^f = <5r7i^^ + 24-^^15-" (4.57) 

= -{Sz^ + 2Sl-d%) . (4.58) 
A more detailed expression for the superpotential is 

= (pA^rii"^]^ - (5mi"^]) ^'^ - m^/'^^A^ - 2C'"t^'^^ ■ (4.59) 
Here the perturbed expressions are used: 

SU^a = U^a-U'ga, (4-60) 

6m'^^ = mf-^t^ (4.61) 

^n)f = n^f-tf- (4.62) 

The perturbation 5z'^ is defined by the the definition (4.10) and by the perturbations (4.61) 
and (4.62). Of course, if a displacement vector is a Killing vector in the background space- 
time then Z^^^ disappears. In the case of GR the conservation law (4.54) goes to the KBL 

conservation law (3.55), if is defined as fc" in (2.20). 

The presented here procedure gives well defined currents and superpotentials in the 
following sense. Without changing the identity (4.7) one can add to the current an arbitrary 
quantity Az"(^) satisfying [Az"(^)]_q; = 0. Analogously, without changing in (4.12) 
the superpotential can be added by A?°'^(^) with the property [A^"'^(^)] ^3 = 0. However, 
the "broken" current and superpotential can be "restored" by the same way because the 
quantities Ai'^(^) and Az"^(,^) are not connected at all with the procedure applied to the 
given Lagrangian in a non-exphcit form. Whereas, in the sense of the procedure, the current 
and the superpotential in Eq. (4.12) are given by the coefficients (4.4) - (4.6) uniquely 
defined by the Lagrangian. This claim develops also the criteria by Szabados [168] who 
suggested to consider a connection of pseudotensors with Lagrangians "as a selection rule 
to choose from the mathematically possible pseudotensors". Thus one can assert that the 
current and superpotential in (4.54) are defined by the unique way in the above sense. Of 
course, the same claim is related to the KBL quantities in (3.55). They uniquely are defined 
by the Lagrangian (3.41) in the sense of the Ncether procedure. 

To construct the Belinfante corrected conserved quantities for the perturbed system 
(4.50) we again turn to subsection 4.1. For the coefficients u^^, m^^ and fi'^^'^ we con- 
struct the identity (4.20) and subtract the same barred identity 



igkO-igkO 



(4.63) 



After substitution the equations (4.52) and (4.53) into the left hand side of (4.63) one gets 
the conservation law: 

i^B){0 = dpi^^){0- (4.64) 
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The current is 

+ (4.65) 
where the Behnfante corrected energy-momentum and Z-term are 

tlB) = -dU^^-Shf,m'p^, + D^5sf,, (4.66) 
^ Sz'^BiO- (4.67) 

The perturbations Sz^^ and Ssg''^a in (4.65) are defined by the definitions (4.22) and (4.19) 
and by the perturbations (4.61) and (4.62). A detailed expression for the superpotential is 

= 2(p/nJ^/l^ + :D,5n;f/1^5pa)r-|^fl^^^^ (4.68) 

In the case of GR the conservation law (4.64) goes to the Belinfante corrected conservation 
law (3.78). We conclude also that in the sense of the united Noether-BeUnfante procedure 
the current and the superpotential in (4.64) are defined in unique way by the Lagrangian 
(4.50). 

Now, rewrite the conservation laws (4.47), (4.54) and (4.64) in the united form: 

iUO=dpiD^iO- (4.69) 
This allows us to construct the conserved charges in generalized form in D-dimensions: 

nO= I d''-'xi%iO=f dSii^i^iO (4.70) 
is ./as 

where S is a spatial {D — 1) hypersurface x° = const and dT, is its [D — 2) dimensional 
boundary. 



4.5 The Einstein-Gauss-Bonnet gravity. 

The mass of the Schwarzschild-anti-de Sitter black hole 

To illustrate the above theoretical results we apply them to the Einstein-Gauss-Bonnet 
(EGB) gravity. The action of the Einstein /^-dimensional theory with a bare cosmologi- 
cal term Aq corrected by the Gauss-Bonnet term (see, for example, [75]) is 

= J {^EGB + jCDm] 

= J d'^x [k-\R - 2Ao) + 7 {Rlupa- ^Rlv + R% + -CDm} (4.71) 

where k, = 2Q.d-2Gd > and 7 > 0; Gd is the D-dimension Newton's constant. The 
equations of motion that follow from (4.71) are 

S^""" = f^"' . (4.72) 
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The metric part is defined as 

+ 27 [RR>''' - iRi'^^R'^P + R^'^prR'"'^ - 2R>'^R'''' 

- \g^''{Rlxpa-^Rla + R^)\} , (4.73) 

In the vacuum case, T^'' = 0, and the equations are 

S^"" = Q. (4.74) 

In the present subsection, as a background we consider the AdS solution, which is a 
solution to the equations (4.74) and is described by the metric: 

D-2 

= -(1 + J)dt^ + (1 + Jy^dr'^ + r'^Y. <labdx''dx'' . (4.75) 

a,b 

The last term describes {D — 2)-dimensional sphere of the radius r, and qab depends on 
coordinates on the sphere only; for the other components ^qo — ~(1 + /) 9ii = 
(1+7)-^ where 

J(r) = . (4.76) 

•'^ ' {D-l){D-2) ^ ^ 

The background Christoffel symbols corresponding (4.75) are 

rJo = ^(1+7)/', ^10 = 7^7/7^^ ^u= 



2(1 + /)' 2(1 + /) 



rib = l^b, rl, = -r{l + f)qab. (4.77) 



The effective cosmological constant (see [75]): 



is the solution of the equation keff~ ^eff^EGB + ^egbM = 0, where 

(D - 2){D - 1) 

is defined only by the Gauss-Bonnet term. Thus Ag// is negative, and the background 
Riemannian, Ricci tensors and curvature scalar are 

R^a.p - 2Aeff p_2)(i^_i) ' ^1^- - ^^-ffD^' ^ - ^^'ffD^ ■ 

(4.80) 

Then equations (4.74) are transformed into 

^M^^ - ha^jR + Ae// V = • (4.81) 
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The Schwarzschild-AdS (S-AdS) solution [215] can be also considered as a solution to 
(4.74): 

D-2 

ds^ = -(1 + f)dt^ + (1 + f)~^dr^ + r^Yl Qabix^dx^ (4.82) 

a,b 

where 



4An 



+ 4k7(D-3)(£>-4)- 



.D-1 



2k7(L>-3)(L'-4) r V * ^EGB 

(4.83) 

For the metrics (4.75) and (4.82) the relation —g = —g = r^~^ det Qab has a place and 
is important for future calculations. The Riemannian, Ricci tensors and curvature scalar 
corresponding to (4.82) are 



-Roioi 

Rahcd 
R 



^ f" , RoaOb — ^^(1 + f)f'Qab , Rlalb — - 



rf 



2(1 + /) 



Qab ■ 



= -r'^fiqacQbd - qadQbc) ; 



2 V" r 
-[f{D-3)+rf']qab; 
,,D - 2 



Rii = 



2(1 + /) 



/" + / 



,D-2 



f" + /' 



^[f{D-3) + rf'] . 



(4.84) 

)■ 

(4.85) 
(4.86) 



Of course, the barred solution (4.82) goes to (4.75), and the barred expressions (4.84) - 
(4.86) go to (4.80). It is evidently, in the case of the solutions (4.75) and (4.82) perturbations 
can be described only by A/ = / — /. In linear approximation it is 



A/ 



4Ao \ 



D-3 



(4.87) 



4.5.1 The field-theoretical prescription 

Here, we turn to the results of subsection 4.3. To concretize and to have a possibility to 
compare with [75], we define the gravitational perturbations from the set h"-, like in (2.77), 

as Kp = gap - 9ap- Thus Ci = Cf^^ = hafsiSCEGs/Sgap) = Kp£ , where £ is 
the barred expression (4.73). Then for the AdS background (4.75) the equations (4.72) can 
be rewritten in the form of the equations (4.35): 



rii^^ — 

^Lq = 



He) 



1 + Wi 



4Ao \ 



Gf^+Jl 



4Ao 

^EGB 



Gi'^ = Kti^^ (4.88) 



where subscripts and (^j^) are related to the Einstein and the Gauss-Bonnet part in 
(4.71), which are with the coefficients "k" and "7", respectively. The left hand side in 
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(4.88) is calculated using (4.36) and 

4 A 



D-2 



2A, 



eff 



(4.89) 



D-2 

The right hand side of (4.88) generalizes the energy-momentum in (2.71). It could be 
presented as 

2 



= -- 



-9 



(4.90) 



where indexes for £^i, and Tf^i, were lowered by g^jii,, and they are thought as depending on 

the sum^^^ + h^p. 

Constructing the charges for the EGB system one has to use the generalized expres- 
sion (4.70). The current could be used. However, looking at (4.48), (4.73) and (4.90) one 
can see that it is very complicated. Evidently that the superpotential expression (4.44) is 
significantly simpler. We set Ci = Cf^'^ in (4.42) and obtain 



-9 
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n 



la — 



=F- 



K 



4Ao _ _ 



4An 



^M^pA ^ ^MP^^A ^ ^\u-p^, _ ^„\-pu _ ^pu-t,\ ^ (gf^^gP^^ _ ^W^^^) . (4.91) 

The substitution of the quantities (4.91) into the expression (4.44) gives the superpotential 
related to the general EGB case: 



4Ao 

^EGB 



X (^['^L>^/i''l'^ - ^^D^hP^'' - ^^''D'^h - h^^^'D'^k^ - ^hD^^(P^^ . (4.92) 

It is expressed through the Abbott-Deser superpotential in the Einstein theory [68, 75], 7^^, 
as X|^'^ = =F-\/l — AAq/ Aegb I ad- Recently, developing the results [75] Deser with co- 
authors [76] have reached the same expression. Paddila [216] reffering to [75] also analyzed 
it. Keeping in mind that different definitions for the metric perturbations /i" can be used we 
find that in fact the superpotential (4.92) belongs to the set 



4Ao \ 



(4.93) 

where the Einstein part defined in D dimensions formally coincides with (3.74) in GR. 

For the S-AdS solution (4.82) considered as perturbations with respect to the AdS 
spacetime (4.75) one has in linear approximation: 
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(4.94) 
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To calculate the total conserved energy of the S-AdS solution we use the formula (4.70) 

where we substitute (Ol)-component of the superpotential (4.92). The last is calculated 
with (4.94) and the time-hke Killing vector A'* = (—1, 0); covariant derivatives are defined 
by (4.77). Finally (4.70) gives: 

^^(Z^2)^ (4.95) 

that is the result of [75], and is the standard accepted result obtained with using the various 
approaches (see [216] - [225] and references therein). 

4.5.2 The canonical Noether mass 

The results of this subsubsection, in fact, repeat the results by Deruelle, Katz and Ogushi 
[182]. However, there is a difference between the approach in [182] and our approach in 
[119]. In [182] conserved quantities are constructed in the way, when the charges for the 
physics system and for the background system are constructed separately. Only after that 
they are compared and differences are interpreted as charges for perturbations. On the other 
hand, in [119] the conserved charges from the starting are constructed in the perturbed 
form. This approach is in the spirit of the present paper and is presented expUcitly here in 
significantly more detail than in [1 19]. 

To construct global conserved quantities in the above prescription we again turn to 
the generalized integral (4.70), only this time with the superpotential (4.59). Then it is 
necessary to calculate the perturbations (4.61) and (4.62). For this one has to use the general 
formulae (4.5) and (4.6) with the metric Lagrangian of the EGB gravity in (4.71). Also the 
expressions in D-dimensions: 

A^. = V^^,-Tl, = \g'=^p{D^gp, + D,gp^-Dpg^,) , (4.96) 

Dpgi^u = gr^^Al^+gruAl^, (4.97) 

i?V = ^pA^a-^aA^^ + A^A?,-A^,A?^ + ;RV (4-98) 
are useful. For calculation of the superpotential (4.59) we need in antisymmetric part of 
(4.61) only. Thus after prolonged calculations we have 

= ^{2A["y]'' + 4"A^]5''^} 

+ 7^ {2i?,''^["A^l - GRr^'^^^PAlp + GRP^'^f^KAlp - 4RP^''A% + 2RPA^^y^^ 
+ ARf-5l^Af!l - 6ii^["/l^A^,ff,, - 124"/-A^, - S^^^^^^A;, + 2RAty^P 



+ 



RA^^X^g"^ + 3i?A;/["/l^ - 2g,pg-^''D,Rf^P - d^^gf^POpR] . (4.99) 
The expression for (4.6) is more simple: 

+ lV^{-2R„^'^^^^ - AR^g^f^ + 4i?^"/)^ + R (g'^^d^ - g^^'^5^^) }(4.100) 
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For calculation of the superpotential (4.59) we need in the next part of (4.100)^: 

- (E)'^a + (GB)K 



E4"/]A + + 4^A[a^^l + ^-/JA^j . (4,101) 



An 

Perturbations of (4.99) and (4.101) are obtained following the recommendation given in 
(4.61) and (4.62): 

= 5[^E)^^f^]+S[^GB)m^f\, (4.102) 
-5 4"/^' = 5[iE)f>^f^^] + ^[iGB)n'^f^^\- (4.103) 

We do not present expUcit expressions for these perturbations because they are evident due 
to (4.99) and (4.101). We note that only the two last terms in (4.99) contribute into n^gf^ 
because ATq,^ = 0. 

Recall that in the canonical prescription one has to define a divergence in the Lagrangian 
(4.50). We exactly follow the recommendation in [182]. In our notations it is the divergence 
of d?" = h^^'^Alp. In the GR case this leads to the choice of KBL, that is to k'^ in (3.40) 
that leads to variation of the action under Dirichlet boundary conditions. A more detail 
discussion on a choice of a divirgence one can found in [182]. Thus, we define 



+ 7^=5 (-2i?<.°^^ + 4i?°/^ - AR^g'^^ + 5^g'''^ - S^g^^) A^^ . (4.104) 

Substitution of (4.102), (4.103) and (4.104) into (4.59) gives the canonical Noether super- 
potential = i^^E ^{^GB ^^'^ perturbations in EGB gravity. Its Einstein part formally 
coincides with the KBL superpotential (3.54): 

^(c)E = -ia'^'^Dpe^ + r^'^A^C - dH^]) - 2^["41 . (4.105) 

only one has to keep in mind that it is presented in D-dimensions [182]. 

To calculate the mass of the S-AdS black hole we use the formulae (4.75) - (4.87) for 
background and perturbed systems. We use the component J^^-^ in the general formula 

(4.70) under the requirement r ^ 00. For this it is enough to calculate Z^^-^ in linear 
approximation with respect to the perturbation A/ in (4.87). Next, for calculating the mass 
we again need in the KilUng vector A'* = {—1, 0}. In linear approximation the symbols 
(4.96) are 

^00 = -Yr^fliD - 5) , A?o = -YrY^D - 1) , 

Aii = l^(Z?-l), A]^ = -rAfq,,. (4.106) 



'in [119] an analogous to (4.100) expression has a missprint. 
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In calculations the simple relations (A/)' = —{D — 3)A//r, (/)' = 2//r are used. Then 
in linear approximation the Einstein part (4.105) gives 



i(c)E = -^^fiD-2). (4.107) 

The contribution from the Gauss-Bonnet part of (4.59) we write out in part. Thus in Unear 
approximation 



-'^[(GB)rnf IJA'^ = 2|_^A//(L> - 2){D - 3) [3 + 2{D - 5)] . (4.108) 



The contribution from (4.103) is 



iA^i^A<5[(GB)^f - |<5[(GB)nf ^'I^aA'^ = -^^A//(i^ - 2)(i5 - 3) . (4.109) 



At last, the contribution from the Gauss-Bonnet part of (4.104) is 



-2A[°d35 = 2^A//(D - 2){D - 3) . (4.110) 

Summing (4.108) - (4.1 10) we obtain the Gauss-Bonnet part of the superpotential (4.59) in 
Unear approximation 



i±Wi 



^EGB 



(4.111) 



The definitions (4.76) - (4.79) were used. Summing (4.107) and (4.111) we obtain the full 
01-component of the EGB canonical superpotential necessary for calculating the mass of 
the S-AdS black hole: 



^S = i^(7j (^-2)- (4-112) 

The definition (4.87) was used. Substitution of (4.112) into (4.70) gives again the standard 
result (4.95). 

4,5.3 The Belinfante corrected mass 



Comparing the Belinfante corrected superpotential (4.68) and the canonical one (4.59) 
we find that for constructing the first we need (together with (4.101)) in the expression 
nYx'^g^^^gpa instead of (4.99) and (4.104). Thus in EGB gravity 

.rp[a /3]A _ 
^gX 9 9 pa — 

+ 4i?^5"[" + ^g'^R^x + R (2'^l5'''" - } . (4.113) 

As is seen, it significantly simplify the expression (4.68) with respect to (4.59). Perturbation 
of (4.113) is again obtained following the recommendation given in (4.62): 

5[n;f /I'V] = S[(E)nY^''9'%a] + SkGB)nfy'^%a] ■ (4-114) 
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Combining (4.103) and (4.114) in the definition (4.68) we obtain the Belinfante corrected 
superpotential X|^^ = i"B)E ■^{b)gb ™ EBG gravity. Its Einstein part 

iii)E = I {i'^"D,f^ + i^'^DpW^P - D^P^e) (4.115) 

formally coincides with the superpotential (3.65), only one has to keep in mind that it is 
presented in L)-dimensions. 

To calculate the mass of the S-AdS black hole we use again the formulae (4.75) - (4.87) 
for background and perturbed systems. We use the component i^B) general formula 
(4.70) under the requirement r ^ oo and for A'* = {—1, 0}. We calculate a linear approx- 
imation of X°g-j in A/. The (01) -component of the Einstein part (4.115) gives the same 
result as in (4.107). The contribution from the Gauss-Bonnet part of (4.68) we again write 
out in part. Thus in linear approximation 

2A-D,<5[(Gs)n7 Vffp.] = ^^AfJ{D - 2){D - 3) [3 + 2{D - 4)] . (4.116) 

Summing (4.109) and (4.116) we obtain the Gauss-Bonnet part of the (Ol)-component 
of the superpotential (4.68) in linear approximation that is exactly (4.111). Finally, thus, 
adding the Einstein part we again approach the standard result (4.95). 



4.6 Discussion 

The results of subsections 4.2 and 4.3, first, generalize the corresponding results in GR 
presented in sections 2 and 3, second, they generahze the Deser and Tekin approach [75] 
in quadratic theories that is particulary presented in the framework of EGB gravity in sub- 
section 4.5.1. Thus really we answer the criticism in [216] related to the Deser and Tekin 
approach. Indeed, in subsections 4.2 and 4.3 we have extended the approach in [75] to 
perfectly arbitrary backgrounds, which can be without symmetries at all, and arbitrary dis- 
placement vectors (no necessity in Killing vectors) can be used. 

The results of subsection 4.4 generalize the corresponding results in GR presented in 
section 3. Also, our canonical Noether expressions generalize the Deruelle, Katz and Ogushi 
results [182] in EGB gravity, which in fact are presented in subsection 4.5.2, but only in a 
preferred here perturbed form. The BeUnfante corrected expressions both in a general form 
in subsection 4.4 and in appUcation to EGB gravity in subsection 4.5.3 are quite new. 

A choice for more suitable and successful expressions for conserved quantities is not so 
simple problem even in the framework of GR. Therefore, only we discuss general propos- 
als for the generalized superpotentials (4.44), (4.59) and (4.68) from a general viewpoint. 
However, a relation to GR, and to EGB gravity (subsection 4.5), of course, could be useful. 

Let us turn to the superpotentials of the family (4.44) with the substitution of (4.42). 
They differ from each other starting from the second order for different I'^^s^ , like the super- 
potentials (3.74) in GR. The choice h"- l^" = g^" - '^'^'^ in GR was based on the two 
points (see the end of subsection 3.5). First, only for this choice the superpotential from 
(3.74) gives the standard energy-momentum and its flux at null infinity for the Bondi-Sachs 
solution. Second, only the superpotential related to Z'^^ from the family (3.74) coincides 
with the Belinfante correction of the KBL superpotential. We do not know a test model for 
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a choice of perturbations /i" in the family (4.44) in EBG gravity. For example, an arbitrary 
choice h"- for the solution (4.82) leads to the same result (4.95) if we use a correspondent 
superpotential from the family (4.93) instead of (4.92). Indeed, for this calculations only 
the linear order is crucial, but it is the same for all of them. Considering the Belinfante 
corrected expression (4.68) and the family (4.44) we see that, in general, they cannot be 
identified. Indeed, the superpotentials (4.44) are linear in l'^^^ for an arbitrary metric theory, 
whereas (4.68) does not. However many modified -dimensional gravities, as a rule, are 
the Einstein theories with corrections (quadratic and others, like string and M-theory cor- 
rections). Among such theories the Lanczos-Lovelock theory (in particulary, EGB gravity) 
plays an important role. Then, comparing the Einstein parts only (in the EGB gravity they 
are (4.1 15) and the one in (4.93)) we again can prefer from the set I'^^y Of course, it is 
not an absolute choice, but could be a recommendation. 

Comparing (4.44), (4.59) and (4.68) we remark that the first is in the framework of the 
field-theoretical approach, where perturbations are examined expUcitly. At the same time, 
the expressions (4.59) and (4.68) are presented rather as bimetric ones. Moreover, the ex- 
pression (4.44), being a linear one, is significantly simpler. Next, the superpotential (4.59), 
as a canonical quantity, essentially depends on a choice of a divergence in the Lagrangian. 
But for an arbitrary metric theory, as we know, there is no a crucial principle for the def- 
inition of such a divergence; there are only reasonable recommendations (see [182] and 
subsection 4.5). On the other hand, expressions (4.44) and (4.68) do not depend on diver- 
gences at all, that could be an advantage. Comparing (4.59) and (4.68) we can also note 
that the Belinfante corrected expressions are simpler significantly. This is demonstrated 
clearly by the corresponding expression in EGB gravity and by the calculation of mass of 
the S-AdS black hole in subsection 4.5.3. 

The appUcations in subsection 4.5 related to calculation of mass of the S-AdS black 
hole in EGB gravity are, in fact, the test of the presented here approach. Recalling the 
papers [216] - [225] where the S-AdS mass was calculated we note that all of them are in 
agreement giving the acceptable result, although they use different approaches. In the most 
of the papers the definitions are based on the first law of the black hole thermodynamics 
and the connection with the surface terms of Hamiltonian dynamics. The surface terms are 
defined following the recommendation by Regge and Teitelboim [190]. On the other hand, 
considering the AdS background as an arena for perturbations one can connect conserved 
charges with AdS symmetries expressed by KilUng vectors and define them in the quite 
classical form. Such approach was used in [182] and in just the present paper. In [224] the 
charges are associated to the diffeomorphism symmetries of Lovelock gravities in any odd 
dimensions. In [225] asymptotic symmetries of AdS spacetime were used with employing 
the conformal completion technique. 

Sometimes it is important to interpret some special situations with the use of developing 
approaches. In the works [76, 216, 221, 226, 227] such a case, when 

Af;GB = 4Ao, (4.117) 

is remarked. Indeed, this case looks as degenerated one. The reason is that the condition 
(4.1 17) leads to a zero coefficient at the hnear approximation of the EGB gravity equations 

around the AdS background. Applying the generalized KBL approach in EGB gravity [182] 
Deruelle and Morisava [226, 227] have found that mass and angular momentum expressions 
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for the Kerr-AdS solution in EGB gravity have the same coefficient and have to be defined 
only as zero. In the framework of the field-theoretical formalism we approach the analogous 
conclusion. Turn to the equations (4.88). Under the condition (4.117) their linear left hand 
side disappears: the Einstein part is canceled by the Gauss-Bonnet part. Consequently the 
total energy-momentum for perturbations at the right hand side disappears also. This means 
that the matter energy-momentum is compensated by the energy-momentum of the metric 
perturbations. In the vacuum case the energy-momentum of the metric perturbations is 
equal to zero. Thus the conserved integrals have to be equal to zero, like the charges in 
[226, 227]. This situation on the AdS background is similar to the definition of zero energy 
of the S-dS black hole in a so-called Nariari space in [222]. The last is interpreted as a 
real ground state whose energy is lower than energy of the pure dS space. To escape the 
difficulties with the situation (4.117) Paddila [216] using his approach suggests to use a 
modified background without all the symmetries, unUke AdS one. All three the approaches 
presented in the present paper have this possibility also. 

In future we plan to check the formulae of the subsection 4.5 applying them to the 
Kerr-AdS solution in EGB gravity keeping in mind, firstly, like in [226, 227], that asymp- 
totically in linear approximation it has to go to the solution in the Einstein D-dimensional 
theory [228, 229, 230], secondly, checking the just now appeared exact solutions [231, 232]. 
Moreover, the expressions related both to the field-theoretical and to the canonical Noether 
prescriptions already were used and discussed for examination of higher D Kerr-AdS space- 
times [76, 226, 227, 230] including EGB gravity. Planning to check the formulae in EGB 
gravity of subsections 4.5.1 and 4.5.2 for the Kerr-AdS solution we anticipate that the re- 
sults will be identical with the ones in [76, 226, 227, 230]. AppHcations of the Belinfante 
corrected expressions of subsection 4.5.3 for calculating mass and angular momentum for 
the Kerr-AdS black hole in EGB gravity have to be testable. 

Gravitational theories in D dimensions are very intensively developed. New solutions 
with interesting properties appear following one by one (see, for example, [233, 234]). To 
understand and describe their properties, which could be quite unusual, many of new solu- 
tions need in calculating conserved charges. Therefore all of known and new possibilities 
to define and calculate conserved quantities, if they are non-contradictive and satisfy all the 
acceptable tests, could be interesting for future studies. 
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